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ABSTRACT 

The  results  of  a  one-year  research  effort  addressing  some  fundamental  scientific  issues 
relevant  to  particle  diagnostics  in  optically-thick  sprays  are  presented.  The  objective  of  this 
research  effort  involved  development  and  application  of  an  expeiiroental  apparatus  for  studying 
scattering  in  optically  thick  me^a.  The  important  technical  contributions  of  this  project 
included  1)  development  of  a  fluidized  bed  approach  for  creating  controlled  optically  thick 
media,  and  2)  characterization  of  performance  of  multi-element  detectors  and  light-valve  arrays. 
With  respect  to  the  former,  we  have  developed  a  binary-particle-phase  fluidized  bed  concept 
combing  with  refractive-index  matching  that  will  allow  the  bed  to  be  fluid-dynamically  very 
dense  (with  interpardcle  spacings  on  the  order  of  two  diameters  as  required  for  stable  operation) 
but  optically  less  thick  (with  interpardcle  spacings  greater  than  five  diamet^  where 
independent,  as  opposed  to  dependent,  muldple  scattering  is  in  effect).  Condidons  for  stable 
operadon  of  the  binary-particle-phase  fluidi^  bed  have  been  identified.  Regarding  detector 
and  light-valve  arrays,  Ae  experimental  results  indicated  diat  edge  effects  are  important  and  that 
current-generation  Faraday-effect  light  valve  arrays  do  not  have  the  performance  specifications 
required  for  use  in  multi-angle  interrogation  schemes  for  diagnostics  in  optically  thick  sprays. 

I.  INTRODUCTION 

Particle  and  droplet  size  distributions,  being  parameters  of  fundamental  importance,  should 
be  priority  measurement  objectives  for  intelligent  sensors  in  next-generation  propulsion  systems. 
Howev^  the  potential  appUcation  of  laser  scattering  systems  as  optical  sensors  introduces  some 
severe  .  equirements  on  the  measurement  techniques.  In  particular,  the  presence  of  multiple 
scattering  in  the  near-field  regions  of  sprays  significantly  complicate  die  application  of  optical 
diagnostics.  This  research  project  involved  experiment^  work  related  to  measurements  under 
such  conditions.  In  particular,  a  fluidized  bed  system  with  the  potential  for  providing  a  stable 
medium  with  controlled  and  continuously  variable  optical  depth  for  a  constant  size-distribution 
was  proposed  and  studied.  The  present  work  was  placated  on  a  preceding  project  as  discussed 
byHirleman  [1,2]. 

The  long-term  objective  of  this  research  is  to  understand  the  process  of  multiple  scattering 
as  it  might  occur  in  future  applications  of  intelligent  particle  size  distribution  sensors  and 
^velop  inversion  schemes  which  can  operate  in  such  an  environment  A  minimum  requirement 
is  that  the  presence  of  multiple  scattering  be  diagnosed  by  the  instrument  to  ensme  that 
erroneous  data  not  be  used  in  control  algorithms,  a  situation  which  could  result  in  a  catastrophic 
failure.  Our  goal  was  to  si^ass  that  and  develop  efficient  and  robust  algorithms  which  can 
actually  extract  useful  particle  size  information  from  Fraunhofer  diffraction  measurements  in 
multiple  scattering  environments. 

The  laser  diffraction  particle  sizing  technique  comprising  the  framework  for  this  research 
has  been  described  in  detail  in  the  literature;  for  current  references  see  [3].  In  the  sections 
below,  the  results  of  this  project  are  discuss^  This  final  report  is  intended  as  an  overview  and 
synthesis  of  the  research  results  as  opposed  to  an  exhaustive  restatement  of  what  has  appeared  in 
previous  progress  reports,  annual  reports,  and  published  technical  papers.  Publications  resulting 
from  this  work  are  incorporated  by  way  of  references  from  the  technical  section.  Finally,  a 
summary  of  personnel  involved  in  the  project  is  presented. 

n.  FLUIDIZED  BED  SYSTEM 

The  creation  of  a  stable  and  well-characterized  two-phase  system  such  as  a  liquid  droplet 
spray  (in  a  gas)  or  a  dispersion  of  solid  or  liquid  particles  in  a  liquid  is  very  difficult  and  has 
escaped  all  efforts  of  standardization.  Probably  th&  most  straightforward  approach  to  obtain  a 
known  calibration  particle  system  for  optical  instruments  is  to  disperse  glass  or  polystyrene 
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calibration  microspheres  available  as  certified  reference  materials  in  a  liquid  using  stirred-  or 
flow-cells.  Our  early  efforts  in  experimental  verification  utilized  such  an  apparatus  but  questions 
concerning  uniform  (spatial)  dispersion,  centrifugal  separation  (by  size)  of  the  particles  in  broad 
size  distributions,  and  the  performance  of  the  system  at  very  high  particle  loadings  (Le.  at  or  near 
the  dependent  scattering  limit)  remained  in  question. 

To  solve  some  of  the  problems  associated  with  cells  such  as  segregation  and  deposition  we 
designed  and  constructed  a  fluidized  bed  system  [4]  as  shown  schematically  in  Fig.  1  and  in  a 
photograph  in  Hg.  2.  In  theory,  a  fluidized  bed  should  provide  a  uniform  dispersion  of  the 
particles  and  therefore  a  very  simple  means  by  which  to  control  the  particle  concentration  and,  in 
turn,  the  optical  depA.  The  fluidued  bed  consists  of  the  particles  of  interest  placed  in  the  test 
section  between  two  microsieves.  A  fluidizing  liquid  is  pumped  in  a  closed  loop  around  the 
system  which  includes  a  filter.  When  the  fluid  passes  through  the  bed  the  resulting  drag  force  on 
the  particles  will,  if  the  velocity  is  higher  than  some  critical  value,  overcome  gravity  and  lift  the 
particles.  If  the  system  is  operating  properly,  for  a  given  fluid  velocity  the  bed  of  solid  particles 
expand  until  the  overall  buoyancy  and  drag  forces  are  balanced.  The  bed  should  expand  and 
contract  as  the  flow  rate  is  varied,  and  the  concentration  (#/cc)  of  particles  will  just  depend  on 
the  height  of  the  expanded  bed  and  the  total  number  of  particles  originally  plac^  in  the 
fluidizing  section. 

The  fluidizing  section  in  the  bed  developed  here  and  shown  in  Figs.  1  and  2  is  7.62  x  2.54 
cm  (3  X  1  inch)  inner  dimensions  with  a  height  of  30.48  cm  (12  inch).  Optical  access  is  provided 
by  two  anti-reflection-coated  quartz  windows  providing  a  2.54  cm  (1  inch)  diameter  clearance 
near  the  bottom  of  the  bed  and  a  7.62  cm  (3  inch)  path  length  for  the  laser  beam.  (The  capacity 
for  similar  windows  on  opposite  sides  of  the  shorter  dimension  is  available).  A  plexiglass 
window  miming  the  full  height  of  the  fluidizing  section  provides  optical  access  to  allow 
measurement  of  the  bed  height  for  the  independent  determination  of  particle  concentration. 

With  the  two  optical  paths  the  system  can  theoretically  allow  a  10:1  variation  in  optical  depth 
while  nrnntaining  the  same  particle  size  distribution  by  merely  adjusting  the  flow  rate.  The 
absolute  values  of  the  optic^  depths  are  determined  by  the  number  (mass)  of  particles  placed  in 
the  fluidizing  section  of  the  bed  when  it  is  assembled. 

We  have  successfully  stabilized  the  operation  of  the  fluidized  bed  under  several 
conditions.  A  photograph  of  the  system  in  operation  is  given  in  Fig.  3.  However,  the  stable 
operating  conthtions  we  have  found  have  been  for  void  fractions  (defined  as  liquid  volume 
divided  by  the  total  volume  of  the  bed)  too  low  for  the  multiple  scattering  models  relevant  here. 
In  particitiar,  90%  void  fraction  is  the  highest  value  for  which  we  have  achieved  stable  operation. 
At  higher  values  of  void  fraction,  the  bed  underwent  the  well-known  "spouting"  phenomenon 
and  the  necessary  uniform  spatial  distribution  of  the  particle  phase  was  lost.  This  level  of  void 
fraction  (90%)  corresponds  to  an  inteiparticle  spacing  of  approximately  2.7  diameters,  a  value 
which  falls  within  the  dependent  scattering  regime  of  Fig.  4.  While  dependent  scattering  is  of 
general  interest,  the  objective  of  this  project,  based  on  multiple  scattering  models  developed  in 
Ae  previous  work  that  apply  only  in  the  independent  scattering  regime,  therefore  requires  higher 
void  fractions. 

To  overcome  the  problem  of  insufficient  void  fraction,  we  propose  a  binary-particle- 
phase  fluidized  bed.  This  idea  is  original,  and  no  reference  to  a  similar  concept  has  been  found 
to  date.  The  solid  phase  of  the  bed  ^1  be  comprised  of  two  discrete  modes,  one  which  we 
designate  as  the  fl  .idizing  particle  medium  and  the  other  as  the  scattering  medium.  These  are 
differentiated  via  the  use  of  a  refractive-index-matching  fluid  which  renders  the  fluidizing 
particle  medium  invisible.  The  scattering  particles  are  of  a  difrerent  refractive  index  and  are 
therefore  visible  and  scatter  light  passing  through  the  bed.  The  scattering  phase  must  have  the 
same  fluid  dynamic  properties  as  the  fluidizing  particle  phase  so  as  to  not  be  segregated  by  the 
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Figure  2.  Fhowgraph  of  the  ASU  fluidized  bed  system. 


Figure  3.  Photograph  of  the  test-section  of  the  ASU  fluidized  bed  system  in  operation.  The  fluid 
was  water,  and  the  particle  phase  was  made  up  of  glass  beads  in  the  74  -  105  diameter  range 
obtained  fiom  Cataphote  Corp.  The  operation  was  stable  at  this  condition,  with  a  bed  height  of 
3.8  cm  (1.5  inch). 


optical  depth 


Figure  4.  Optical  property  map  for  a  fluidized  bed.  Curves  of  constant  physical  path  length  1  are 
shown  vs.  inter-particle  spacing  and  optical  depth  for  100  lim  spheres.  The  regimes  of 
dependent  and  independent  scattering  are  shown  as  well. 
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fluidized  bed.  In  this  way  the  void  fractions  of  the  two  solid  particle  phases  are  independently 
variable,  thereby  allowing  control  of  the  optical  depth  of  the  medium  independent  of  the  total 
void  fraction.  In  this  way  the  bed  can  be  operated  in  a  fluid-dynamically  stable  regime  while 
optically  in  the  independent  scattering  regime  as  is  required  for  this  work. 

The  design  problem  basically  then  involves  selection  of  the  properties  of  the  two  particle 
phases  and  the  fluid.  One  critical  constraint  is  that  the  fluidizing  particle  phase  must  be  optically 
homogeneous,  i.e.  no  inclusions  or  other  inhomogeneities  inside  Ae  particles  can  be  tolerated. 
We  have  had  difficulty  obtaining  particles  of  acceptable  quality.  In  order  to  achieve  refractive- 
index  matching,  it  is  easier  to  obtain  fluids  to  match  the  refractive  index  of  a  given  particle  rather 
than  vice-versa.  To  check  for  a  refractive  index  match  we  developed  a  procedure  based  on  the 
Beck's  Line  method  that  uses  refraction  phenomena  near  the  edge  of  a  particle  to  qualitatively 
ascertain  the  ratio  of  the  refractive  indices  of  the  particles  and  the  fluid.  Figures  5  and  6  show 
photographs  taken  through  a  microscope  (40X  objective)  of  glass  beads  in  the  74-105  pm 
diameter  range  obtained  from  Cataphote  Corporation  immersed  in  Cargille  fluids  of  refractive 
index  1.51  and  1.52  respectively.  In  Fig.  5,  the  outline  of  the  particles  are  quite  visible 
indicating  that  the  specified  nominal  refractive  index  of  the  particles  of  1.51  is  off.  Also  clearly 
visible  in  Fig.  5  are  a  number  of  impurities  and  inclusions  within  the  spherical  particles  that 
remain  visible.  These  inclusions  would  scatter  light  in  an  experiment  and  therefore  render  these 
particles  unacceptable  for  the  binary-particle-phase  bed.  Figure  6  is  a  similar  photograph  where 
the  particles  have  been  immersed  in  a  fluid  of  refractive  index  1.52  where  excellent  matching  has 
occurred.  The  outlines  of  the  particles  are  effectively  invisible,  though  the  inclusions  within  the 
particle  remain  visible. 

Due  to  the  rather  poor  optical  quality  of  these  particles  shown  in  Figs.  5  and  6,  it  was 
necessary  to  look  elsewhere  for  particles.  Some  example  of  spheres  of  better  quality  are  shown 
in  Figs  7  and  8.  Probably  the  best  quality  particles  we  have  tested  are  the  high-index  spheres 
from  Duke  Scientific  Coi^.  shown  in  Fig.  8.  The  problem  with  these  is  the  high  refractive  index 
of  1.91,  since  the  available  refractive  index  fluids  in  that  range  are  highly  toxic  (Arsenic 
compounds)  and  certainly  very  difficult  to  use.  Pyrex  particles  obtained  ftom  Mo-sci  near  the 
end  of  the  project  have  reasonably  good  quality  and  are  shown  in  Fig.  9. 

Our  proposed  design  of  the  binary-particle-phase  fluidized  bed  would  involve  Mo-sci 
Pyrex  or  Yttrium-Alumina-Silica  beads  as  the  fluidizing  particle  medium  matched  with  Cargille 
fluids  and  the  high  density  glass  microspheres  from  Duke  Scientific  as  the  scattering  medium. 
For  stable  operation  in  the  ASU  fluidized  bed  for  optical  depths  of  1  - 10  we  propose  600  pm 
diameter  particles  for  the  fluidizing  particle  and  420  pm  diameter  for  the  scattering  medium. 

As  of  the  conclusion  of  this  project  we  have  been  unable  to  demonstrate  the  binary- 
particle-phase  fluidized  bed  concept.  The  particles  and  fluids  required  for  our  proposed  design 
are  estimated  (based  on  retail  pricing  at  small  quantity  prices)  of  the  order  of  $10,000. 
Demonstration  of  our  binary-particle-phase  fluidized  b^  concept  will  be  the  subject  of  future 
work. 

m.  CHARACTERIZATION  OF  DETECTOR  AND  LIGHT  VALVE  ARRAYS 

In  the  previous  project  [1]  we  suggested  a  possible  scheme  for  a  general  solution  to 
Fraunhofer  diffraction  particle  size  measurements  in  multiple  scattering  envirorunents.  The 
experimental  system  to  allow  multi-angle  interrogation  of  the  particle  field  is  shown  in  Fig.  10. 
The  key  difference  between  Fig.  10  and  conventional  single  scattering  Fraunhofer  instruments  is 
the  presence  of  a  programmable  mask  in  the  ftont  focal  plane  of  the  transmitter  lens.  The 
programmable  mask  has  aimular  ting  apertures  which  can  be  individually  switched  on 
(transmitting)  or  off  (absorbing  or  opaque).  A  ring  of  light  in  the  front  focal  plane  of  an  ideal 
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Figure  5.  Photograph  of  particles  in  matching  fluid  taken  through  an  optical  microscope  with  a 
40X  objective  and  a  red  filter.  The  particles  (of  diameter  nominally  one-third  of  the  field-of- 
view)  are  glass  beads  in  the  diameter  range  74-105  pm  obtained  from  Cataphote  Corp.  The  glass 
beads  have  a  specified  nominal  refractive  index  of  1.51,  and  they  are  dispersed  in  a  matching 
fluid  of  refractive  index  1.51  obtained  from  Cargille. 


Fig^  6.  Photograph  of  particles  in  matching  fluid  taken  through  an  optical  microscope  with  a 
40X  objective  and  a  red  filter.  The  particles  (of  diameter  nominally  one-third  of  the  field-of- 
view)  are  glass  beads  in  the  diameter  range  74-105  pm  obtained  from  Cataphote  Corp.  The  glass 
beads  have  a  specified  nominal  refractive  index  of  1.51,  and  they  are  dispersed  in  a  matching 
fluid  of  refiractive  index  1.52  obtained  from  Cargille.  The  refractive  index  matching  here  is  quite 
good,  and  basically  only  the  inhomogeneities  inside  the  particles  are  visible. 
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Figure  7.  Photograph  of  particles  in  matching  fluid  taken  through  an  optical  microscope  with  a 
40X  objective  and  a  red  filter.  The  particles  (of  diameter  nominally  one-half  of  the  field-of- 
view)  are  Yttrium-Alumina-Silica  beads  of  nominally  100  pm  diameter  obtained  from  Mo-sci 
Co.  The  YAS  beads  have  a  specified  nominal  refractive  index  of  1.69,  and  are  dispersed  in  a 
matching  fluid  of  refractive  index  1.69  obtained  from  Cargille. 


Figure  8.  Photograph  of  particles  in  matching  fluid  taken  through  an  optical  microscope  with  a 
40X  objective  and  a  red  filter.  The  particles  (of  diameter  nominally  one-sixth  of  the  field-of- 
view)  are  high  density  glass  microspheres  obtained  from  Duke  Scientific.  The  glass 
microspheres  have  a  specified  nominal  refractive  index  of  1.91,  and  are  dispersed  in  a  fluid  of 
refractive  index  1.54  obtained  from  Cargille. 
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Figure  9.  Photograph  of  particles  in  matching  fluid  taken  through  an  optical  microscope  with  a 
40X  objective  and  a  red  filter.  The  particles  (of  diameter  nominally  one-fourth  of  the  field-of- 
view)  are  Pyrex  spheres  obtained  from  Mo-sci  Co.  The  spheres  have  a  specified  nominal 
refractive  index  of  1.49,  and  are  dispersed  in  a  fluid  of  refractive  index  1.54  obtained  from 
Cargille. 


Figure  10.  Schematic  of  optical  system  for  multi-angle  interrogation  of  optically  thick  sprays. 
The  programmable  mask  on  the  input  side  transmits  cylindrical  shells  of  light  which  are 
converted  to  hollow  cones  of  light  by  a  lens.  /  matched  programmable  mask  at  the  transform 
plane  allows  me^urement  of  the  redistribution  matrix.  The  insert  shows  the  actual  geometry  of 
the  34-element  ring-aperture  light  valve  designed  and  fabricated  according  to  the  optimal  scaling 
law  developed  by  Hirleman  [5]  in  the  preceding  project. 


transmitter  lens  produces  a  hollow  cone  of  light  of  constant  angle  6  passing  through  the  spray  or 
particle  sample  volume.  By  switching  open  die  various  rings,  a  sequence  of  hollow  cones  of 
probe  radiation  is  created.  The  fraction  of  the  incident  energy  in  the  cones  which  is  not  scattered 
by  the  medium  is  redirected  by  the  ansform  lens  to  a  ring  on  the  detection  plane  which  matches 
Ae  ring  in  the  programmable  mask  (assuming  the  focal  lengths  of  the  transmitting  lens  and  the 
transform  lens  are  equal).  Light  which  is  scattered  by  particles  in  the  spray  leaves  at  some  angle 
different  thw  the  cone  angle  and  ends  up  at  another  radial  position  (Le.  a  different  detector  ring) 
on  Ae  detector  plane. 

The  method  proposed  for  generating  the  cones  of  light  necessary  for  implementing  the 
multi-angle  interrogation  scheme  propos^  by  Hirleman  [2]  for  deconvolution  of  the  multiple 
scattoing  effect  involves  the  use  of  Faraday-effect  light-valve-arrays.  Three  custom  devices 
have  been  fabricated  based  on  the  optimal  scaling  geometry  of  Hirleman  [S].  Initial  efforts  at 
experimental  verification  are  report^  by  Kenney  and  Hirleman  [4].  The  success  of  the  method 
hinges  on  die  performance  of  detector  and  light  valve  arrays.  The  contrast  ratio  must  be  such 
that  there  is  minimum  cross-talk  between  elements.  Our  efforts  at  utilizing  these  devices  have 
basically  been  unsuccessful.  The  problem  appears  to  involve  the  performance  of  die  individual 
elements.  It  was  necessary  to  perform  a  careful  characterization  of  the  performance  of  these 
devices  in  order  to  isolate  and  quantify  the  problem. 

The  experiments  used  in  the  portion  of  the  work  were  similar  to  those  reported  by  Kenney 
and  Hirleman  [6,7,8].  A  laser  beam  was  focussed  onto  the  light  valve  array  and  a  photodetector 
was  placed  in  line  with  the  beam  behind  die  light  valve  to  collect  the  transmitted  light  The 
results  from  a  scan  across  two  detector  elements  of  the  light-valve  array  are  shown  in  Hg.  1 1  for 
the  light  valves  electronically  switched  on  (open)  and  then  off  (closed).  The  contrast  ratio 
observed  for  this  experiment  is  about  6,  si^uficantly  less  than  that  observed  in  work  supported 
during  the  previous  projert  by  Hirleman  and  DcUenback  [9],  It  appears  that  the  performance  of 
the  devices  has  degreed  over  the  past  year,  and  is  now  at  a  stage  where  the  cones  of  light 
generated  by  the  device  are  not  adequate  for  the  application.  Significant  advances  in  the 
technology  will  be  required  before  Faraday-effect  devices  will  be  usable  in  the  multi-angle 
interrogation  scheme  proposed  here. 

Results  of  some  additional  characterization  experiments  on  the  arrays  are  shown  in  Fig.  12. 
Here  the  focussed  laser  beam  (nominally  30  pm  diameter)  was  scanned  across  the  entire  array  by 
moving  the  array  in  2.5  pm  st^s.  Also  shown  in  tire  lowest  Fig.  12  are  similar  results  indicating 
the  resolution  of  the  experiment  obtained  by  scarming  beam  across  a  knife  edge.  It  is  seen  that 
the  edge  defirution  on  tiie  light  valve  array  is  considerably  better  tiian  diat  for  photodiode  array 
as  reported  by  Kenney  and  Hirleman  [7,8].  The  uniformity  across  the  element,  as  indicated  by 
the  fluctuations  in  transmitted  energy  at  the  tops  of  the  waveforms  in  Figs.  1 1  and  12  is  similar 
to  that  we  measured  for  photodiode  arrays  [7,8]. 

In  sunmiary,  the  performance  of  the  light  valve  arrays  fabricated  for  this  project  was  not 
adequate  for  incorporation  into  multi-angle  interrogation  schemes  for  diagnostics  in  optically 
thick  sprays.  At  present,  until  technology  advances  improve  the  performance  of  the  light  valves 
significantly,  other  means  must  be  used  to  generate  the  interrogation  beams.  Future  work  will 
use  axicons  or  leflexicons  to  this  end 
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Figure  12.  Transmitted  optical  energy  vs.  position  as  a  focussed  HeNe  laser  beam  was  scanned 
across  the  ASU  light  valve  array  ^m  the  center.  A  knife-edge  scan  of  the  beam  (see  data 
indicated  by  A  on  the  left  side  of  the  lower  figure)  indicated  a  beam  diameter  (l/e^  intensity 
points)  of  30  pm .  The  edge  definition  of  the  light  valves  edges  is  much  better  than  30  jim.  The 
transmission  indicated  through  the  inner  rings  (uppo*  figure)  does  not  saturate  since  the  rings  are 
thinn^  than  the  laser  beam. 
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Vm,  SUMMARY  OF  IMPORTANT  TECHNICAL  CONTRIBUTIONS 

1.  Development  of  the  binary  particle-phase  fluidized  bed  concept  This  idea,  combined  with 
refractive-index  matching,  will  allow  the  bed  to  be  fluid-dynamically  very  dense  (with 
interparticle  spacings  on  the  order  of  two  diameters  as  required  for  stable  operation)  but 
optically  less  thick  (with  interparticle  spacings  greater  than  five  diameters  where 
independent  as  opposed  to  de^ndent  multiple  scattering  is  in  effect). 

2.  Experimental  characterization  of  the  performance  (edge  effects,  uniformity,  and  contrast 
ratio)  of  programmable  optical  shutter  arrays  for  use  in  multi-angle  interrogation  of 
optically  thick  media.  References  [6,7,  and  8]. 


From  Proceedings  of  the  2nd  International  Congress  on  Optical  Particle 
C  Sizing,  March  5-8,  1990,  Tempe,  AZ. 

T  *  A  GENERAL  SOLUTION  TO  INVERSE  FRAUNHOFER  DIFFRACTION 

PARTICLE  SIZING  IN  MUL'nPLE  SCATFERING  ENVIRONMENTS:  EXPERIMENT 


Steve  B.  Kenney  and  E.  Dan  Hirieman 
Mechanical  and  Aerospace  Engineering  Department 
Arizocu  State  University 
Tempe.  AZ  85287-6106 

One  of  llie  more  difficiilt  remaining  challenges  in  the  field  of  optical  particle  sizing  is  the  characterization  of  particle  size  distributions 
of  optically  iltick,  multiple  scattering  media.  In  a  previous  paper  [1],  a  theoretical  approacli  to  solving  \i»  forward  prtAlcm  (i.e. 
calculation  of  the  scattering  signature  given  the  particle  size  distribution)  for  the  case  of  near-forward  scattering  by  p^des  large 
compared  to  the  wavelength  in  optically  tliick  media  was  presented.  A  subsequent  paper  included  in  tliis  Congress  12]  is  concerned 
witli  die  corre-spoiiding  inverse  problem,  i.e.  tlrat  of  estimating  the  size  distribution  function  of  an  ensemble  of  large  particles  using 
madiematical  inversion  of  the  near-forward  optical  scattering  properties  of  llie  optically  thick  medium.  In  this  paper  we  discuss  the 
experimental  approach  wliicfa  will  be  used  to  verify  the  dKoretical  developments. 

Tlie  inverse  .scattering  problem  for  single  scattering  requires  at  least  n  measurements  (wliich  is  equivalent  to  specifying  n  equatioas)  to 
solve  for  n  unknowas,  where  die  unknowns  arc  tire  quantity  of  particles  in  n  size  classes.  The  n  measurements  are  die  amounts  of 
light  scattered  into  various  scattering  angles.  Conventional  Fraunhofer  diffraction  systems  use  a  collimated  (incident  angle  of  zero) 
interrogation  laser  beam,  but  incident  light  at  any  near-forward  angle  could  in  fact  be  used.  The  optical  system  of  Fig.  1  generates 
itxrident  (conical)  beams  of  various  angles,  and  for  each  incident  cone  n  scattering  measurements  (i.e.  at  n  detectors)  could  be  inade. 
Titus,  the  system  in  Fig.  I  could  theoretically  provide  enough  measurements  to  support  n  independent  solutions  to  the  n  by  n  single 
scattering  system.  The  optical  system  in  Fig.  1  system  uses  transmission  mode  spatial  liglit  modulators  as  discussed  by  DcUenback 
and  Hirieman  t3]. 

In  optically  thick  media,  the  inverse  problem  is  significantly  complicated  because  the  li^t  which  evenmally  reaches  the  detectors  has. 
in  geiKral,  undergone  more  than  one  scattering  event.  In  odicr  words,  the  single  scattering  signature  (which  can  be  modeled  relatively 
easily)  is  perturbed  or  altered  by  tlie  additional  (ntultiple)  scattering  events  which  occur  after  the  light  leaves  the  first  scattering  event 
on  its  way  to  a  detector.  Additional  unknowns  are  thereby  iiHroduced,  specifically  variables  which  determine  or  predict  the  expeaed 
fate  of  light  which  would  have  reached  each  of  the  various  detectors  if  not  for  the  intervening  thick  medium.  Now  since  light 
originally  traveling  at  a  .scattering  angle  corresponding  to  say  the  jth  deteaor  could  be  rescattered  into  n  other  detectors,  one  can 
define  n^  unknowns  necessary  to  model  the  n-angle  scattering  beliavior  of  an  optically  thick  medium.  In  contrast,  the  single 
sc.'ittering  medium  required  only  n  equations  and  n  unknowns.  OearW  then  addidonaJ  measurements  are  needed  to  characterize  a 
multiple  scattering  mediutn,  and  the  system  in  Fig.  1  can  provide  independent  measurements  (each  of  n  incident  light  cones 
.scattered  to  n  detectors).  After  the  n^  measurements  are  taken  via  Fig.  1,  an  or  by  n^  system  of  equations  must  effectively  be  solved 
to  detennine  the  particle  size  distribution. 

It  has  been  slwwn  [2)  (hat  the  inverse  problem  for  multiple  scattering,  in  the  case  of  near-forward  scattering,  can  involve  the  on-line 
measurement  of  the  multiple  scattering  redistribution  matrix  H,„  defined  by; 

(1) 

where  Sq  is  a  vector  defining  the  irradi.tnce  of  incident  light  in  n  angles,  and  S,„  Ls  a  vector  defining  tlie  angular  distribution  of 
irradiance  of  the  (multiple-scattered)  light  exiting  tlie  medium.  Tlie  scalteriiig  redistribution  matrix  H,p  such  tliat  tlie  matrix  element 
in  the  ith  row  and  the  jth  column  H^fi.j)  is  tlie  gain  or  efficiency  witli  which  optical  energy  incident  in  the  jth  angle  or  direction  is 
redistributed  by  tlie  medium  into  the  itli  direction.  Now  we  also  know  [2]: 

1  exp(-b)  •  expfa^b  •  H)  (2) 

where  af  and  b  are  the  forwaid-scatteting  albedo  and  the  optical  depth  respectively  and  H  is  a  n  by  n  redistribution  matrix  valid  for 
single  scattering.  Now  a  conventional  diffraction  system  under  multiple  scattering  conditions  will  measure  a  signature  proportional  to 
tlie  first  column  of  H,„,  and  only  for  small  b  are  H  and  H„,  equivalent.  But  H  can  be  determined  from  measured  H  using: 

H  =  exp(b)/(a(b)  •  ln(  )  (3) 

Tlie  n^  temis  of  II^,  are  measured  using  Eq.  (1)  by  illuminating  tlie  medium  using  n  different  Sq  vectors  (i.e.  n  different  hollow  cones 
of  incident  illumination)  and  measuring  tlie  scattering  signature  on  n  different  detectors  for  each  incidence  case.  Then  H  is  obtained 
from  Eq.  (3)  whidi  in  effect  amounts  to  "reacliing  into  the  medium”  and  detemiining  the  scattering  signature  tliat  is  present  after  just 
one  scattering  event  lias  taken  place.  Once  the  single  scattering  properties  of  the  medium  are  known,  the  inverse  problem  can  be 
solved  using  an  array  of  methods  discussed  in  the  literature;  see  Koo  [4]  for  references. 
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Experinienial  verificalioo  of  the  algorithm  for  inverse  scattering  under  multiple  scattering  conditions  requites  an  e^tmental 
apparatiB  where  Uk  optical  density  and  particle  size  distribution  of  the  medium  can  be  carefully  comroUed.  In  both  the  nui^cal  and 
experimental  studies  performed  here  we  use  NBS  standard  reference  material  (SRM)  1003A  as  the  particle  size  distributioa  In  the 
fusi  expcriiiienis  we  used  a  optical  cell  with  a  magnetic  sdning  system,  and  d^ised  the  SRM  1003a  particles  in  water  in  die  cell 
Tlic  tlieorctical  optical  density  was  proportional  to  particle  concentration  assuming  that  the  panicles  were  uniformly  distributed  in  the 
cell.  Particle  concentration  in  tlie  liquid  was  vari^  to  produce  panicle  fields  of  different  optical  thicknesses.  Stirred  cells  do  have 
potential  problems  in  tliat  tlie  centrifuge  effect  can  cause  both  size  and  spatial  segregation  of  the  particles. 

The  generation  of  hollow  cones  of  Illumination  is  performed  using  a  spatial  light  modulator  operated  in  the  transmission  mode.  The 
custom  device  has  ihiny-two  semicircular  ring-shaped  apenures  which  can  be  individually  addressed.  Tlie  ring  geometries  foUow  die 
optimal  scaling  law  discussed  by  Hirleman  [5].  Activating  any  ring  element  causes  tliat  element  to  rotate  the  polarizatitm  vector  of 
tlie  transmitted  light  relative  to  that  of  light  passing  through  the  unactivated  elements.  A  polarizing  element  beliind  the  ring  window 
array  then  preferentially  blocks  light  which  has  not  been  rotated.  This  in  effect  produces  conical  sliells  of  light  passing  through  the 
tratismission  ina.sk,  which  in  turn  produce  conical  shells  of  liglit  at  the  back  focal  plane  of  the  lens.  Cycling  tlirough  the  ring  elements 
produces  a  sequence  of  cones  of  light  of  different  included  angles.  A  matched  device  at  tlie  back  focal  plane  of  the  collection  lens 
allows  measurement  of  die  n^  elements  of  H. 

Particle  Flow  System 

A  fluidized  bed  was  designed  to  create  a  sample  volume  in  which  the  particle  concentration  could  be  easily  varied.  A  sdiematic  of 
tlie  experimental  apparatus  is  shown  in  Fig.  2.  IruUally,  the  particles  to  be  sized  lie  on  a  stairdess  steel  mesh  at  the  bottom  of  the 
fluidizing  tube.  If  a  fluid  in  laminar  flow  is  passed  upwards  through  the  tube  and  past  the  particles  a  critical  velocity  will  be  reached 
wlien  die  viscous  frictional  and  drag  forces  equal  the  weight  of  the  particles  in  the  fluid  stream.  As  the  flow  rate  is  increased  the  bed 
expands.  For  beds  of  liquid-solid  composition  die  bed  expands  unifonnly  and  the  concentradoii  of  solid  particles  is  evenly  distributed 
throughout  tlie  bed  L6J.  Therefore,  varying  optical  depths  can  be  produced  simply  by  changing  the  flow  rale  of  the  fluid. 

A  0.2nim  fllter  is  placed  in  the  line  to  keep  unwanted  foreign  material  from  entering  the  fluidizing  tube.  To  achieve  a  uniform  flow 
across  the  crosis-seclional  area  a  packed  bed  of  glass  beads  (6.S4niin  diameter)  was  inserted  before  the  fluidizing  tube  entrance. 
Optical  windows  with  anti-reflection  coatings  are  mounted  near  the  bottom  of  the  fluidized  tube  on  opposite  sides.  The  irKident  laser 
light  is  passed  tlirough  tlie  window  and  tlie  light  scattered  by  the  suspended  particles  is  focused  and  collected  on  the  other  side.  A 
rectangular  window  running  tlie  length  of  tlie  fluidizing  tube  is  used  to  detcmiine  the  height  of  tlie  bed  from  wliicii  the  particle 
couceuuation  con  be  calculated.  To  ensure  unifonii  flow  all  windows  are  mounted  flush  to  the  inside  surface. 

Results 

Experimental  results  obtained  witli  the  multiple  scattering  cell  apparatus  discussed  above  are  shown  in  Fig.  3.  Here  the  measured  data 
for  two  columns  of  tlie  matrix  are  plotted  against  the  calculated  (tlieoretical)  signatures  coupled  with  the  successive  order 
multiple  scattering  model.  The  data  .show  reasonable  agreement,  and,  based  on  the  theoretical  inversion  studies  in  the  presence  of 
noise  1 1 1,  it  apfiears  that  recoiLstructions  for  reasonable  numbers  of  degrees  of  freedom  will  be  possible.  Experiments  to  characterize 
the  full  H„|  matrix  are  mxlerway. 
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Pig.  I.  SchcmMic  of  i  nexi-geneiwian  Icwr  (tiffnaiim  panidc  tiring  imnnHiiem  which  allow*  multi-angle  imenogalion  and  multi-angle 
tcalleritig  meaauiemeias.  TIk  progfanimaUe  mask  creaiea  atauilar  ring  of  light  wlucti  tit  convened  into  cone  of  liglil  of  varying  angles 
to  iiluniinate  tiie  panicle  ncki. 


Fig.  2.  Sdiemiic  of  Ouiiiized  bed  system. 


ng.  3.  Mcaanrtd  and  predicted  data  for  two  eolomns  of  the  muhiple  actutering  redistribution  mairia  H  ^  The  data  are  for  MBS  SRM 
IU03A  for  a  bner  diffraction  syaiem  with  an  KSI  Inc.  ring  detector  and  a  ttwitfonn  tens  with  f  ■  123  mm.  Tbe  j  values  coneapood  to  ring 
detector  numben  10.0162  and  0.(M34  rad  for  ring*  14  and  22  respectively). 
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ABSlRACr 

Photodiode  anays  used  in  laser  difihaction  particle  sizing  insuuments  most  be  calMwaled  to  aocoont  for  delector-to-delector 
varialioos  in  sensitivity.  We  have  calibialed  an  Insilec  EPCS-P  ring  detector  Qdanar-diffiised,  p  on  n  silicoo  photodiode 
array)  by  scanning  a  focussed  laser  beam  across  foe  detector  saifooe.  A  deconvolnlioo  of  foe  known  intensity  dhiriibation  of 
foe  laser  beam  from  foe  measured  signal  resulted  in  detector  response  as  a  frmclion  of  poailion.  Detector  reqwose  was 
spptozimaldy  constant  over  foe  region  of  foe  ideal  active  detector  and  h  decreased  exponentially  in  foe  region  beyond  foe 
ideal  detector  boundaty.  A  difinskm  length  constant  of  SOpaa  gave  foe  best  fit  to  foe  measured  data.  Iheoieiical  piedktiaos 
of  caUbration  factors  baaed  on  measured  detector  leyonse  agreed  reasonably  wdl  with  Malvern  and  Insilec  cafibration 
fociofs  obtained  from  foe  iracfitional  unifrmn  light  illmninrtioo  mefood.  This  mdicates  that  edge  effects  in  (fifierent  ting 
detectors  are  similar. 


Many  optical  processing  systems  require  an  amy  of  pholodetectors  to  measure  photon  fttx  at  the  focal  |dane.  We  defoie  a 
detector  as  a  device  that  converts  electromagnrtictadiilion  into  an  dectrical  signal.  The  magnitude  of  the  dectrkal  signal  is 
dfaecdy  mlated  to  the  photon  flux  incident  on  the  device.  Photodelectors  have  been  commercially  available  for  some  time  in 
1>D  and  2*D  arrays  of  various  elemeni  shapes.  Ritnre  improvements  in  system  performance,  from  a  detector  standpoht,  will 
include  larger  arrays,  incieaacd  speed  and  dymanic  range,  and  higher  resoluikiL^  Higher  lesohaioo  transiaies  into  smaller 
elemems,  and  foeir  assocUled  problems  such  as  (1)  crosstrik,  foe  generation  of  signals  at  elements  other  than  the  one  being 
iHuminated  and  (2)  sensitivity,  or  foe  unequal  response  of  elemems  to  foe  same  electromagnetic  flux.  On  a  smaller  scale, 
dlffercnoes  in  etetnem  to  efetnem  sensitivity  result  ditecriy  from  foe  spatial  notmnifr>tm  reqxmse  over  each  demem  area.^ 
This  lespcnse  ncmunifrMinity  is  due  to  oonditiaiis  such  as  imperfectioos  in  foe  crystalline  latliw,  variatioiv  in  doping  density, 
dominam  modes  of  efectron-hde  transport  (diflusion  or  dtiftX  and  surfine  defects  and  contamination. 

To  accoum  for  nommifiMniilies,  detector  siptds  are  often  corrected  by  mnll^ying  foe  signal  by  a  predetennioed  constant 
These  conectioo  fisctois  are  determined  empiiicdly  by  exposhig  foe  detectors  to  one  or  more  nnifoim  itradiaooe  levels.^ 
This  procedure  compensates  for  rdadve  pixd  to  pixd  sensitivity  differences,  but  does  nm  yield  any  infbrmalioo  on  qwlial, 
imra-pheel  reaponsivfty  variations.  Oeneraily,  foe  most  dramatic  variarion  in  qwtid  reqwosivity  oocurs  near  foe  e^  or 
boundary  of  the  detector.  As  dements  become  smaller,  foe  tesponsivhy  near  foe  efetnem  edges  has  a  mote  pronounced 
effect  on  foe  overall  sensitivity  of  foe  demern.  The  "edge  effect”,  or  responeivity  near  foe  demem  edge,  is  the  topic  ttf  ftus 
paper,  partfeidaily  m  to  how  te  edge  eflect  impacts  log-scded,  antwiar,  semi-ciicular  shaped  detector  arrays. 

An  example  of  an  qtpUcalioo  of  phmofoode  turays  is  foeir  use  in  particle  sizmg  instruments^  based  on  near  frwwatd 
scaneting  signaimes.  Scattering  signatures  are  cdlecied  by  a  detector  array  and  are  processed  by  a  mafoemadcd  invetsioo 
scheme  fo  detenrine  the  size  disttibmioa  Embedded  in  foe  soinlioo  is  an  assumed  detector  array  lerponse.  Consequendy, 
one  foe  Ihnfts  on  foe  accnacy  of  the  aiae  disitibolioo  sdutioo  is  proper  esfonatian  of  te  detector  arr^  sensitivity. ' 

In  foe  inetiuinem  described  drove  a  system  matrix  is  uauaily  assembled  assuming  a  step  detector  response  defined  by  the 
detector  geometry.  An  improvemem  in  the  mdtix  accuracy,  however,  would  include  a  mote  teaharic  detector  reyrose. 
Through  «qKtimem  we  have  measured  foe  response  of  a  dmector  at  various  qrotid  poims  with  the  aid  of  a  focussed  laser 
beam,  ftom  foe  exprrimcntd  ttota  we  have  defined  a  mathematical  modd  that  accuniidy  predicts  detector  response  through 
foe  boundary  region.  Abo,  we  show  how  caKbratioo  factors  are  highly  dependem  on  er^e  terrorise  and  that  bom  a 
foeoieticd  model  of  foe  edge  response  dwracteristice  it  fa  possible  to  predict  calibratioofirotota. 
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A  photodiode  detector  is  manufiKtuied  from  a  tfain  siiicoa  wafer.  One  ride  of  the  wafer  is  doped  with  atoms  from  the  tfiird 
column  of  the  periodic  table  such  as  boron,  while  the  other  side  is  doped  with  atoms  from  the  fifth  ccdumn,  such  as 
phosphorus.  Silicon  doped  with  Boron  is  minus  an  electron  to  balance  the  charge  in  the  crystalline  structure  and  the  material 
is  called  p-type  silicon.  This  gives  the  material  a  slight  positive  charge,  hence  the  symbol  p  is  used,  likewise,  silicon  doped 
with  Phosphorus  is  termed  n-type  silicon  becaitse  of  the  extra  electron  in  the  lattice,  giving  it  a  slight  negative  diarge. 

Between  the  p  and  tHype  silicon  regions  exists  the  p-n  jurxtion.  When  the  two  types  of  silicon  are  brought  together  the  extra 
electrons  from  the  n-type  material  diffuse  to  the  p-type  to  fill  up  the  holes  near  the  Boron  atoms.  The  junction  has  a  width  of 
otdy  a  few  atoms  and  a  small  potential  is  establish^  due  to  the  extra  positive  diarge  at  the  n  boundary  and  the  negative 
duurge  at  the  p  boundary. 


Silicon  rioms  readily  absorb  photon  energy.  If  the  absorbed  energy  is  greater  than  l.leV,  an  electron-hole  pair  will  be 
oeated.  Pairs  creried  in  die  depletion  region  or  within  a  diffiisioo  length  constant  will  be  separated  by  the  electric  field, 
leading  to  current  flow  in  the  external  circuit.  For  a  p-on-n  semicondudor,  holes  will  colled  at  the  from  ohmic  contad, 
while  electrons  will  cdled  at  the  back  ohmic  comad.  The  resulting  curtent  is  approximately  proportional  to  the  tdal 
illumination  and  because  of  this  relationship  the  photodiode  can  be  used  as  a  light  measuring  device.  (For  additional 
information  on  the  physics  of  phmodiode  detedors,  see  Sze).^ 


3.  SIZING  INSTRl 
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r  MATHEMATICAL  BACKGROUND 


Scattering  rneasureniems  in  laser  diffraction  instruments  are  tyincally  made  with  annular  ring  sluqied  detedors  which  cover  a 
finite  range  of  scattering  angles  as  determined  by  the  detedor  qiertures.  It  is  conveniem  to  assign  a  particular  scattering 
angle  B,  to  tepresem  the  range  of  scattering  angles  in  the  a^rerture.  The  detection  process  can  be  represented  as: 


iw(0i)=/ we(e,ei)i(e)de  d) 

0 

where  the  weighting  function  W0(6,Oj)  describes  the  relative  responsivity  variations  across  the  detedor  and  I^B})  is 
tepresemative  of  the  sigrud  obtain^  from  the  ith  discrete  detedor  (the  subscript  w  indicates  a  dependence  on  the  weighting 
function).  The  azimuthal,  or  9  dependence  on  the  responsivity  variations  are  significantly  less  important  than  the  0  effects, 
and  have  been  neglected.  i(6)  is  imensity  (W/sr)  scattered  at  near-forward  scattering  angles  6.  The  weighting  function 
wq  is  what  is  to  be  detentrined  experittMmally.  It  is  r^o  noted  that  the  particle  distribution  should  also  be  discretized  in  a 
manner  similar  to  that  used  for  the  scattering  imensity. '  In  that  case  we  obtain  a  system  of  m^equations  in  m^fa  is  the  size 
parameter)  unknowns  where  m^  is  the  number  of  dtsmte  detedors  and  m^  is  the  number  of  discrete  size  classes.  The  linear 
system  is  written  as: 


I  =  K-N  (2) 

In  Eq.  (2)  the  mg  elements  of  the  vedor  I  are  I^6j)  as  given  by  Eq.  (1);  the  m^  elements  of  N  contain  the  bth  partial 
moments  of  the  number  of  p^cles  in  the  size  clam;  and  K  is  the  instrumem  or  system  matrix  whereby  elemem  K;: 
represents  the  diffraction  contribution  of  a  unit  measure  of  partides  in  the  jth  size  clam  onto  the  ith  ddedor.  The  elements 
of  K  ate  given  by: 


Kjj-/  /kb(a,0i)w„(a,ctj)wg(e,ei)dode  (3) 

00 

where  w^fcuoL)  is  a  weighting  fundion  for  the  jth  size  class  and  k{^<z,0|)  is  a  general  scattering  function  which  gives  the 
scittlermg  contribution  of  a  unit  quantity  of  particles  of  size  a  imo  angle  0.  Now  the  sdution  or  measured  partide  size 
distribittion  indicated  by  N  in  Eq.  (2)  can  in  thrary  be  obtruned  by  inverting  the  matrix  K. 
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The  eiqienmeiilal  appaialw  consisted  of  a  laser  difCraction  particle  sizing  instrument  (Hg.  1).  A  laser  beam  is  q>atiaDy 
filtered,  eq>anded  and  collimated  to  a  few  millimeter  1/e^  diameter.  Particles  within  the  collimated  beam  scatter  light  which 
is  collected  by  a  receiving  lens  and  refiacted  onto  the  detectw  plane.  A  commmctaUy  available  Insitec  detector  and 
associated  electronics  were  used  in  the  erqreriment.  The  layout  of  the  Insitec  detector  is  weO-suhed  for  examining  detector 
response,  given  its  odd  numbered  rings  on  one  half  of  the  detector  and  even  numbered  rirtgs  on  the  other  (Fig.  2).  This 
design  has  large  noo-reqronsive  sections  between  adjacent  rings,  which  provides  good  isolation  between  aryacent  detecttns. 

The  Insitec  detector  was  mourned  in  a  carKilevered  position  on  a  prommmable  x-y  translation  stage  which  moved  in  a  plane 
normal  to  the  beam  axis.  A  HeNe  laser  beam  was  focussed  to  a  Ife^  spot  25pm  in  diameter  and  the  qxH  was  positiotied  in 
the  small  hole  in  the  center  of  the  detector.  (The  hole  allows  the  strong  unscattered  light  at  the  center  of  the  defector  to  pass 
through,  thereby  reducing  stray  reflections  and  inner  ring  cross-talk),  (hoarse  centeririg  was  performed  by  moving  the 
detector  until  the  light  reflecting  off  the  region  surrounding  the  hole  disappeared  (indicating  the  li|^  was  passing  through  the 
hole).  The  beam  power  was  adjusted  so  that  when  the  laser  was  positioned  in  the  middle  of  an  active  defectm  the  resultant 
signal  was  just  below  saturation  level.  Several  tuns  were  made  by  scanning  across  the  entire  detector  on  a  line  passing 
through  the  center  in  one  ditectiao  and  then  rotating  the  defector  90  degrees  and  repeating  the  scan  in  the  other  directioa 
Signals  were  recorded  at  Ipm  increments  and  the  entire  experimental  apparatus  was  stelded  from  stray  rrxrm  light. 

Before  scmning  the  defector,  the  tings  were  examined  under  a  nucrosoqre.  Three  distinct  regions  were  noticed  which  ate 
shown  in  Hg.  3  and  ate  labeled  as  follows:  the  ideal  responsive  or  active  region  (Rgctive^*  ^  transition  region  (gap),  and  the 
noo-active  tegioo  (nuaked).  The  p-n  junction  within  the  active  region  is  also  shown,  though  not  to  scale.  The  non-active 
region  between  each  defector  consists  of  a  thin  alumitaim  film  which  has  been  deposited  to  mask  the  area  between  detectors. 
It  is  assumed  that  the  metal  helps  block  radialioo  so  that  masked  regions  are  not  reqKnsive.  Though  Hg.  3  shows  just  one 
detector,  all  detectots  were  examined  and  found  to  be  similar. 

To  measure  the  gap  width,  the  defector  was  placed  on  a  two-axis  translation  stage  that  was  coupled  to  LED  displays  showing 
relative  stage  position.  The  tesolittion  of  die  display  was  Ipm.  A  camera  connected  to  a  19-iDch  TV  screen  was  mounted  on 
the  microscope  which  allowed  for  easy  visual  impection  of  the  object.  A  small  piece  of  tape  was  placed  on  the  screen  to 
mark  the  starting  and  ending  points  for  various  features.  An  inner  gqi  edge  was  visually  lined  up  with  an  imaginaiy  tangent 
line  on  die  edge  that  would  be  parallel  to  the  tape.  The  sta^  were  moved  until  the  outer  gap  edge  was  reached  and  the 
difference  between  stage  positions  was  recorded.  In  this  marxier  the  width  of  each  gap,  (both  iimer  and  outer  per  detector) 
was  determined. 


5.  RESULTS  AW  DISCUSSION 

The  average  gap  width  was  approximately  8pm  (results  from  all  defector  rings  can  be  found  in  Table  1).  The  largest 
measured  gap  width  was  10pm  while  the  smallest  gap  width  was  6pm,  though  most  widths  were  8pm.  The  table  begins  with 
defector  3  because  the  first  four  inner  rings  were  damaged  when  a  hole  was  laser  drilled  through  the  defector  center.  Also 
included  in  Table  1  is  the  visual  width  of  the  Insitec  detectors  (from  microscope). 

Hg.  4  displays  the  individual  rii%  signal  vs.  radial  distanoe  for  the  first  five  inner  odd  tings.  The  rectangular  boxes  represent 
the  specified  width  of  the  detectots  and  therefore  the  ideal  top  bat  response  distribution,  while  the  curves  represent  the 
measured  signal.  The  measured  signal  is  a  convolution  of  the  defector  ^atial  reqwnse  and  the  Gaussian  intensity  profile  of 
the  laser  beam.  Since  the  measured  curves  are  centered  around  the  themtical  points,  it  was  concluded  that  the  laser  beam 
was  approximafely  centered  about  the  detector,  which  also  confirms  the  accuracy  of  the  stage  motion.  The  maximum  signal 
is  dose  to  8  volts,  however  this  value  is  not  quite  reached  on  rings  4  and  5  because  the  width  of  these  rings  is  less  than  twice 
the  beam  width,  (where  twice  the  beam  widfo  represents  more  than  99%  of  the  total  beam  energy).  Also,  overlap  between 
rings  occurs  becmne  of  the  beam  whMi  and  edge  effects  resulting  from  electron  bole  pair  and  photon  (fiffiisioa  Note  how 
overlap  decreases  with  increasing  ring  size,  which  is  due  to  the  larger  regions  of  non-responsive  area  between  rings.  The 
carves  ^ipear  smooth  except  for  the  sfeuration  region  where  sroaU  deviations  in  signal  are  noticed.  It  should  be  mentioned 
that  1000  signate  are  avera^  at  each  positioa 
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nf.3.  C>(Mi witioa of aalfMilec photodiode deieooraiiowiiififae ideal tetpoosive 
icfioa  (R-aclive),  the  tnneitioa  iqpoo  (gtp),  and  the  noo-active  legioa  (nuaked). 


nf.4.  Plolofdeiectoriipialva.iadiaidiMaaoe  Cor  five  inoer  odd  tinp.  Ibeieciaa^ 
icpeeacat  the  specified  dnector  dfaseniiooe.  The  oarwa  aee  te  Mpdfancnddag  iron  a 
2SfMi  dbOMier  iBKr  beam  of  wavdeagdi  0.6328|ini.  MeaMveaMnia  were  made  al 
l|im  incTraMott  and  dwconeipaodiiii  signal  at  each  poeiiioaii  an  avetafe  of  1000  dabi 
pohtt. 
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Insiiec  ideal  detector  width  and  gap  wkMi.  Visual  dinwnsiotia  wete  acquhed  by  viewing  the  detector  under  a  microacope. 
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In  Older  n>  be  able  to  calcidaie  w^  for  nae  in  Eq.  (1)  k  is  necessary  to  model  die  gq>  responae.  This  teymae  can  be 
eadmaied  by  deconvolving  the  measured  signal  from  the  known  Oannian  imeiisity  pro^  of  the  laser  beam.  After  several 
mns  were  taken  of  the  inner  odd  rings  (results  of  lings  1 1,13,13)  an  average  lesponae  of  the  gap  r^on  was  cdculated  by 
fint  deteiminkig  the  average  signal  for  each  ring  at  each  Ipm  posidon  (over  a  lOOpm  dntance)  and  then  superimpoaiiig  the 
profiles  from  the  3  outer  edges  (S  runs  for  each  edge)  to  cdculale  an  overaD  average.  This  avenge  gap  lesponae  curve  is 
shown  in  Pig.  5. 

Inhia)  insight  into  a  possible  model  for  the  gap  response  is  found  by  looking  for  symmetiy  widun  the  curve.  Hg.  6  shows 
two  curves  labeled  upper  and  lower.  The  upper  curve  is  the  first  SOpm  of  the  avenge  edge  reqxmse  owe  normalised  to 
one,  while  the  lower  curve  is  obtained  by  "inverdog"  the  iMtSOpim  of  the  avenge  edge  curve.  The  lower  curve  is  found  by 
subfraeting  1  from  each  lower  signal  and  superimpoaiiig  the  new  curve  onto  the  upper  curve.  It  is  apparent  from  6  that  a 
certain  d^iee  of  symmetry  esials.  This  is  eapedaily  true  fm  the  signal  values  aasochuedwith  a  posidon  between  30|im  and 
S0|im.  For  raifial  posMons  less  than  30|im  diere  is  sli^  disciepancy  between  awes  which  is  most  likely  attributable  to 
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nndom  noiae  or  reflectioiis  from  the  mask  aEEecting  the  ring  output,  although  the  background  signal  tvfaicfa  contains  most  of 
the  noise  has  already  been  subtracted  from  the  total  signal.  There  is  a  residual  signal  obtained  the  beam  is  illnminating 
the  mask  even  though  we  would  expect  no  teqxNise. 

The  symmetry  present  in  the  curve  points  to  a  potential  model  based  on  a  sim|de  function,  most  likely  eiqxxiential.  A  logical 
model  to  preset  the  gap  response  would  include  uniform  response  over  the  idealized  detector  region,  decaying  exponential 
response  over  the  8pm  giq)  width  and  no  re^ranse  in  the  masked  or  non-iesponsive  region.  Photons  iixxdent  on  the  gap 
region  are  absorbed  by  the  silicon  where  they  create  electron-hole  pairs.  The  electron-hole  pairs  begin  diffusing  through  the 
material  until  they  either  reach  the  p-n  junction,  where  there  is  a  possibility  that  they  will  contribute  to  the  overall  current,  or, 
they  lose  their  energy  through  ccdlisioos  and  recombine.  As  far  as  this  model  is  cotmmed  only  one  parameter  will  need  to 
be  determined,  the  "diffusion  length  constant"  governing  the  decaying  exponential. 

However,  it  should  be  mentioned  that  when  convolving  the  Gaussian  beam  with  the  gap  response  function  the  value  of  the 
diffrision  length  constant  has  only  a  slight  impact  on  the  slope  of  the  curve,  thereby  making  it  difficult  to  estinute  the  prr^r 
value  of  the  constant  to  match  the  measured  data.  Hg.  7  confirms  this  observation,  where  three  curves  representing  the  same 
beam  diameter,  but  different  diffusion  length  constants  are  displayed.  As  can  be  imagined  the  gap  response  resulting  from 
difierent  diffusion  constants  is  hi^y  variable,  yet  because  the  beam  width  is  much  greater  than  the  grp  width  the  slopes  of 
the  curves  are  nearly  equal.  If  the  beam  width  was  less  than  the  edge  width  then  the  slopes  of  the  three  curves  would  be 
much  di^rent.  Rirthermore,  if  the  beam  was  reduced  to  a  point  source  then  the  response  would  simply  be  the  exponential 
itself.  However,  even  with  large  beam  widths  a  differerxre  exists  between  the  curves.  The  results  are  sensitive  to  diffusion 
length  in  that  the  curves  expand  with  increasing  diffusion  length  constant. 

Since  the  qiecified  width  of  each  detector  is  known,  a  theoretical  calculation  could  use  diffusion  constants  as  a  means  to 
expand  or  contract  the  curve  to  fit  the  measured  data  and  determine  the  diffusion  length  constant.  First,  the  total  laser  power 
in  the  Gaussian  laser  beam  is  normalized  to  one  and  the  beam  is  moved  in  one  micron  itKrements  over  the  ideal,  gap,  and 
masked  regions.  At  each  position  the  expected  signal  is  calculated  by  convolving  the  reqxmse  function  (assuming  a 
diffusion  length )  with  the  Gaussian  beam  intensity.  (As  can  be  seen  from  Hg.  2  there  are  instances  when  the  width  of  the 
beam  covets  all  three  regions  of  the  detector.)  These  values  are  then  compared  to  the  measured  data  to  find  the  diffrision 
length  constant  that  provides  the  best  fit. 

Experimental  data  from  ring  9  were  compared  to  theoretical  calculations  based  on  distent  diffrisioo  length  constants. 
Assumed  values  in  the  theoretical  calculations  include  a  2SMin  Gaussian  intensity  profile  diameter  laser  beam,  gap  width  of 
8pm,  and  total  detector  width  obtained  from  visual  inaction.  Results  showed  that  a  diffrisioo  length  constant  of  SOpm  gave 
the  best  fit  between  experimental  arxf  theoretical  values.  A  plot  of  ring  9  normalized  sigiul  vs.  distance  is  found  in  Hg.  8, 
and  as  can  be  seen  frpom  the  gnph  the  fit  is  quite  good.  Difteion  lengths  are  fundamental  characteristics  of  semicoodnetor 
devices  and  numy  of  the  techniques  used  to  measure  them  are  reviewed  by  Schroder.^  Our  measured  diffrision  length  of 
SOpm  ±  lOpm  is  consistent  with  typical  values  for  p-n  junction  devices.  As  a  test  for  the  model,  comparisons  were  made 
between  eiqierimental  and  measur^  results  for  many  of  the  other  tings  and  similar  results  were  noted.  Thus,  it  was 
concluded  that  a  simple  expoiKntial  function  iqipeats  to  adequately  predict  the  edge  reqionse  of  the  Inshec  ring  detector. 
The  weighting  function  over  the  width  of  an  edge  can  be  written  as: 

'^eedge  • 

where  p(6)  is  the  local  responsivity,  q,  is  the  maximum  reiqioiisivity  in  the  ideal  active  region  of  the  detector,  xd  is  the 
diffrision  length  constant,  a^  x  is  the  position  measured  from  the  ideal  edge  into  the  g^).  Now  that  the  weighting  function 
wq  has  been  determined  a  new  instrument  matrix  K  can  be  calculated  by  substituting  the  measured  gap  weighting  function  of 
Eq.  (4)  into  Eq.  (3)  in  place  of  the  usually  assumed  step  fiinctioa  To  quantify  the  effect  of  the  detector  weighting  function 
on  die  instrument  matrix,  two  matrices  were  calculated;  one  assuming  a  step  deteaor  reqwnse,  the  other  the  actual  measured 
response  (and  weighting  values  of  0  and  1  in  the  masked  and  idealized  regions,  respectively).  A  plot  of  the  difference  matrix 
(element  by  element)  between  the  two  matrices  is  presented  in  Hg.  9.  The  matrices  were  calculated  on  a  particle  area  basis 
with  rectangular  weighting  functions,  uniform  by  volume  wiihin-class  distributions,  and  a  Friuti^fer  diffraction 
approximation.  The  surface  plot  shows  almost  no  change  on  the  outer  rings,  while  the  di^rence  between  elements  increases 
as  ring  dimensions  decreases.  This  is  expected,  as  the  ratio  of  edge  area  to  visual  detector  area  is  much  greater  for  smaller 
rings  dian  larger  ones. 


SP/f  Vol.  1 480  Sensors  and  Sensor  Integration  119911  /  89 


0  so  100  1B0 

racM  dManoo  ^tm) 

Hg.  7.  Plot  of  theoretical  signal  curves  calculaied  by  con*iriviQg  a  2Siini  beam  diameier 
with  the  lhe<Metical  detector  reqxmse.  Ihe  specified  detector  width  is  55|un  wfaicfa 
coneqwods  to  ring  9  of  the  Inshec  detector.  Calcularions  are  nuKL!  with  different 
diffusion  length  constants  as  indicated.  The  ideal  active,  gap,  and  masked  regions  are 
also  indicated  on  the  figure. 
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Hg.  8.  Comparison  oi  theoretical  and  measured  signal  vs.  relative  distanoe  for  ring  9  ot 
die  Insilec  detector.  The  theoretical  signal  was  calcnlated  baaed  on  a  2S|im  beam 
diameter,  8|im  edge  vridth,  SOpm  dUIusioo  lengtii,  and  SS|im  detector  width. 
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It  has  already  been  mentioned  that  it  is  possible  to  predict  cdibiatioa  ftcton  based  on  g^>  or  transition  region  lesponse  as 
opposed  to  systematic  between-detector  variations  in  the  le^tonsivity.  The  accuracy  of  the  "gap"  calculated  fretois  is  based 
on  the  premise  that  the  major  contribution  to  calibration  ^ois  is  from  gap  leapoose.  A  predicted  caUbiation  firetor  is 
defined  here  as  the  ratio  of  the  ideal  detector  ou^  to  the  actual  detector  ou^t  when  illuminated  with  ligfai  of  uniform 
intensity.  The  ideal  detector  output  is  calculated  by  integrating  die  incident  light  intensity  over  the  ideal,  visual  detector  area 
while  assenting  duti  the  weighting  fiinction  w^  is  constant  and  equal  to  one.  The  actual  detector  output  is  found  in  the  same 
manner  except  the  area  over  winch  the  integration  is  perfonned  irKhides  the  gap  region  and  the  associated  weighting  fiinction 
^Oedge*  In  equation  form  we  have: 


^Tred 


ildA 

/wgldA 


f  "’Sedge  +  f  ‘*^ide«l) 


(5) 


where  Cp,^  is  the  predicted  calibration  factor  and  I  is  the  intensity  (W/m^. 


Traditionally,  calibration  factors  have  been  determined  based  on  illuminating  the  detectors  with  uitifoim  light  intensity.  We 
determined  two  sets  of  calibration  factors  for  the  Inshec  detector.  One  set  of  fiictors  was  calculated  using  Eq.  (5),  while  the 
other  set  was  determined  using  the  uniform  light  method.  Results  of  the  two  methods  are  displayed  in  Rg.  10  along  with 
calibration  factors  from  3  different^*^^’^  ^  Malvern  2600  instruments  determined  empirically  with  the  uniform  light  method. 
Since  the  InsHec  detector  geometry  differs  slightly  from  the  Malvern  detector  geometry,  the  calibration  factors  are  {dotted  vs. 
detector  width.  As  is  seen  from  the  plot  the  curves  are  in  reasonable  agreement.  Noticed  trends  are  that  all  factors  decrease 
in  value  as  detector  area  decreases.  Consequently,  the  inner  detectors  over-reqtond  relative  to  the  outer  detectors. 
Calibration  factors  determined  from  Eq.  (S)  are  limited  by  the  fact  that  experimental  data  were  collected  along  two  lines 
running  through  the  center  of  the  detector  instead  of  scanning  the  entire  detector  surface.  Additionally,  in  Eq.  (5)  it  was 
assumed  that  the  wei^uing  function  equals  a  constant  over  the  ideal  detector  region,  however,  Rg.  4  shows  di^  variation 
in  this  region.  One  other  area  of  uncertainty  which  we  did  rtot  investigate  is  the  ends  of  each  detector.  Though  we  accounted 
for  the  end  regioo  in  the  calculations  the  region  is  different  from  tire  rest  of  the  {tetimeter  regions  in  that  a  masked  area  is  not 
found  8{un  from  the  ideal  detector  boundary.  A  drawback  to  the  uitifotm  li^  method  is  that  linear  detector  response  is 
assumed  with  tM  cross-talk.  If  either  of  the  above  assumptions  is  ix>t  ctmect  then  errors  in  calibnuion  factors  will  {wopagate 
through  the  remaining  frKtors  that  need  to  be  determined.  This  results  from  the  requireurent  that  values  acquired  from 
different  light  levels  most  be  s{tiioed  together.^  Overall,  it  a|)|)ean  that  calibration  factors  calculated  using  Eq.  (5)  comfiare 
reaaondtly  well  with  factors  determined  with  the  uniform  ligix  method,  erqrecially  considering  the  fret  that  hctots  are  for 
different  detectors  from  different  laboratories. 


6.C0MCLUS1QWS 

In  conclusion,  the  transition  region  response  has  been  measured  and  modeled  for  an  Insitec  detector.  A  sim|rie  ejqwnential 
with  a  diffusion  length  constant  of  S0{un  accurately  predicted  the  experimental  data.  The  detector  weighting  frmction  has 
been  calculated  and  incorporated  into  the  instrumeri  matrix  K.  Calibraticm  factors  for  the  Insitec  detector  have  been 
(nedicted  based  on  edge  effects  and  were  compared  with  calibration  fiictors  derived  from  the  uniform  light  source  technique. 
Both  methods  indkaied  that  the  inner  rings  over-respond  relative  to  the  outer,  larger  tings.  It  is  difficult  to  conclude  which 
method  {troduces  the  more  accurate  results.  The  ui^otm  light  method  has  the  disadvantage  that  delectin'  signals  resulting 
from  several  light  levels  must  be  measured  and  the  calibration  curves  s|Hioed  together.  The  gap  res|XNise  model  has  the 
disadvaittage  that  it  is  based  on  2-dimensiooal  results  and  it  must  be  assumed  that  giq>  res|>onse  is  axisymmetric.  Another 
independent  technique  for  calculating  calibration  factors  could  help  verify  whether  the  gap  or  uniform  tight  source  is  more 
accurate.  The  effect  of  using  the  theoretical  calibration  factors  instead  of  the  unifrirm  li^t  source  factors  would,  in  ^neral, 
translate  into  dightly  larger  particle  sizes  because  of  the  higher  signals  on  the  inner  rings. 
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Rg.  9.  Plot  of  the  diffinetioe  matrix  between  based  on  a  stq>  response  detector 

weighting  function  and  K  obtained  with  theoretical  detector  response  weighting  function 
(which  accounts  for  detector  response  in  the  gap  region  between  the  mask  and  idealind 
detector  boundary). 
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Rg.  10.  Rot  of  calibration  fKtms  vs.  detectw  widtii  for  various  instruments  and 
catibration  techniques.  The  numbers  in  the  legend  correqwnd  to  the  following:  1) 
Malvern  2600,  instrumem  #1,  urtifoim  incoherent  illumination.^’^  ^  2)  Malvern  2600, 
instilment  #2,  foctory  supplied  calibration  factors.^  3)  Malvern  2600,  instnunent  #3, 
mtiform  iiKobeient  illomimuioo.l^  4)  Insitec,  uniform  incoherent  ilhminatioa  5) 
Insilec,  dworetical  predictkms  based  on  gap  response  method. 


92  /  SPI£  Voi  1 4S0  Sensors  and  Sensor  Integration  (1991) 


8.  REFERENCES 


1.  J.  N.  Lee  and  A.  D.  Raher,  "Device  Developments  fw  Optical  Infonnation  Processing."  in  Advances  in  Electronics 
and  Electron  Physics,  V.  69,  P.  W.  Hawkes,  Ed.  Academic  Press,  San  Diego,  CA  (1987),  pp.  1 15-173. 

2.  J.  Campos,  A.  Coirons,  and  A.  Pons,  "Response  Unifonnity  of  Silicon  Pliotodiode,"  Aj^Hed  Optics,  V.  27.  pp. 
5154-5156(1989). 

3.  '  L.  G.  Dodge,  "Calibsation  of  the  Malvein  Paiticle  Sizer,"  Applied  Optics,  V.  23,  pp.  2415-2419  (1984). 

4.  E.  D.  Hirieman,  "Particle  Sizing  by  Optical,  Nooimaging  Techniqnes,"  in  liquid  Paiticle  Size  Kfeasurement 
Teciiniqiies.  STP  848,  J.  M.  Tislikoff,  R.  D.  Ingebo,  and  J.  B.  Kennedy,  Eds.  (American  Society  of  Testing 
Materials,  Philadelpliia,  1984),  i^.  35-60. 

5.  E  D.  Hirieman,  "Unceitaiities  in  Matrix  Formulations  of  the  Fraunhofer  Difihaction  Particle  Size  ProMem., 
Conference  of  the  Institute  of  liquid  Atomizatioo  and  Spray  Systems,  (1989),  Irvine,  CA. 

6.  S.  M  Sze,  Physics  ef  Semiconductor  Devices,  John  Wiley  and  Sons,  New  York  (1981). 

7.  E.  D.  Hirieman,  "Optimal  Scaling  for  Fmmhofer  Diffiaction  Particle  Sizing  Instruments,"  Particle 
Characterization.,  V.  4,  pp.  128-133  (1988). 

8.  D.  K.  Schroder,  Semiconductor  Material  and  Device  Characterization,  John  Wiley  and  Sons,  New  York,  (1990). 

9.  G.  S.  Samuelaen,  Univeisity  of  California,  Irvine,  CA  (private  communication). 

10.  E  A.  Hoveruc,  NASA  Lewis  Research  Cerker,  Cleveland,  Ohio  (private  communication). 

11.  L.  G.  Dodge,  Southwest  Research  Institute,  San  Aikonio,  Texas  (private  communicatioa). 


SPI£  Vol.  1 480  Sensors  and  Sensor  Integration  (1991/  /  93 


Rom  Proceaiing  of  the  Fifth  International  Conference  on  Liquid  Atomizarion  and  Sorav  Systems 
NIST  Special  Publication  813,  U.S.  Govcnnent  Printing  Office,  Washington,  DC,  1991. 
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ABSTRACT 

Photodiode  arrays  used  in  laaer  diffraction  particle  sizing  instruments  must  be  calibrated  to  accouitt  for  detecttv* 
toKietector  variations  in  sensitivity.  We  have  calibrated  an  Insitec  EPCS-P  detector  by  scanning  a  small  laser  beam 
across  the  detector  surface.  A  de^volution  of  the  known  intensity  distribution  of  the  laser  beam  from  the 
measured  signal  resulted  in  detector  response  as  a  fuiKtion  of  position.  Detector  response  was  approximately 
constant  over  the  region  of  the  ideal  active  detector  aixl  it  decreased  exponentially  in  the  region  beyond  the  ideal 
detector  boundary.  A  diffusioo  length  constam  of  SOpm  gave  the  best  fit  to  the  measured  data.  Theoretical 
predictions  of  calibration  factois  based  on  measured  detector  response  agreed  reasonably  well  with  Malvern  and 
Insitec  calibration  factors  obtained  from  the  traditional  uniform  light  illumination  method.  This  indicates  that  edge 
effects  in  ring  detectors  made  by  different  sources  are  similar. 

INTRODUCTION 

Particle  sizing  instruments  based  on  near  forward  scattering  signatures  are  commonly  used  to  measure  droidet 
size  distribudons  in  sprays  [1].  Scaaering  signatures  are  colleaed  by  the  instrument  and  are  processed  by  a 
mathemadcal  inversion  scheme  to  determine  the  size  distribudon.  The  equation  which  models  the  scattering 
signature  is  a  classical  Fredholm  integral  equadon.  This  integral  equation  is  often  sqjproximated  by  a  linear  system 
of  discretized  equations  and  the  calculated  size  distributicm  is  only  an  sqrproximation  of  the  actual  distributiotL  The 
cotreqKXiding  coefifideta  matrix  resulting  from  discretization  contains  elements  that  must  be  estimated.  Two  of 
these  components  are  the  distribution  of  partide  size  within  each  size  dass  and  detector  response  characteristics. 

The  ultimate  accuracy  of  the  calculated  size  distribution  is  limited  by  the  accuracy  to  which  these  components  are 
properly  estimated. 

Since  it  is  the  partide  size  distribution  that  is  sought,  it  is  doubtful  that  any  information  on  within  dass  size 
distributions  would  be  known  a  priori.  However,  it  is  possible  to  experimentally  probe  the  detector  teqxmse. 
Literature  references  [2-4]  discwn  calculation  of  the  scattering  matrix  assuming  an  ideal  stq)  response  at  the 
boundary  of  detector  elements.  Any  deviation  in  the  actual  detector  edge  response  characterises  from  this  assumed 
ideal  behavior  would  result  in  a  bias  error  in  the  matrix,  arxl  hence,  also  in  the  calcuiated  size  distributions.  Rirtber, 
all  scattering  matrix  cakulations  repotted  to  date  have  assumed  no  detector-to-deteaor  variSon  in  re^XMisivity,  the 
ideal  situation.  However,  since  intSaboritoty  studies  [5]  highlighted  the  significant  effect  of  between-detector 
respouivity  variSons  on  the  overall  accuracy  of  laser  diffiaction  instruments,  calibration  factors  have  been 
introduced  to  correct  for  these  between-detector  effects.  The  process  of  calibrating  a  photodetector  array  using 
unifimn  ilhimirtatioo  [6]  can  experimentally  correct  for  two  effects,  1)  spatial  variations  in  the  local  re^onsivity 
(ampsAvatts)  dther  within  or  bMween  detectors,  and  2)  edge  effects.  Qearly,  details  of  these  local  effects  ate 
maadeed  by  the  averaging  effect  of  ordfiotm  flood  iUumination.  In  this  piper  we  are  concerned  with  the  details  of 
detector  array  characteristics  on  a  spatial  scale  modi  smaller  than  the  detector  dimensioos.  In  particular,  we  rqwit 
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meMaements  of  local  teqmmivity  both  within  wid  at  the  boondaiiet  of  the  detector  demcnlt.  We 
have  oaed  the  leeultt  to  develop  a  theoretical  model  for  the  edge  effecu  and  we  therefore  ate  able  to  tbeoictically 
predict  calibnaioo  factoia. 


MATHEMATICAL  FORMULATION 

measurements  in  laser  diffraction  in8tn]m«its  are  typically  matte  with  annular  ring  idufted  detecton 
which  cover  a  finite  range  of  scattering  angles  as  deteimined  by  the  detector  q>ertures.  It  is  convenient  to  assign  a 
pfirtiraiiar  scattering  angle  6|  to  represent  the  range  of  scattering  angles  in  the  ith  deteaor  aperture.  The  detection 
process  can  be  lepreseoted  as; 


Iw(0i)“/ we(e.ei)i(0)de  d) 
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where  the  weighting  function  W0(6,6p  describes  die  relative  responsivity  variations  across  the  detector  and  I«j(6^  is 
representative  of  the  signal  obtained  mm  the  ith  discrete  detector  (the  subscript  w  indicates  a  dependence  on  the 
weighting  function).  The  azimuthal,  or  9  dependence  on  the  responsivity  variations  are  significantly  less  important 
than  the  6  effects,  and  have  been  negleaed.  i(6)  is  the  intensity  (W/sr)  diffracted  at  near-forward  scattering  angles 
0.  The  weighting  function  w^  is  what  is  to  be  determiiwd  eiqieriment^y.  It  is  also  noted  that  the  particle 
distribution  should  also  be  discredzed  in  a  manner  similar  to  that  used  for  the  scattering  intensity  [7].  In  that  case 
we  obtain  a  system  of  m9  equations  in  m^  (a  is  the  size  parameter)  unknowns  where  m9  is  the  number  of  discrete 
detectors  and  m^  is  the  number  of  discrete  size  classes.  The  linear  system  is  written  as: 
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(2) 


In  Eq.  (2)  the  m9  elements  of  the  vector  I  are  lw(0j)  as  given  by  Eq.  (1):  the  mg  elements  of  N  contain  the  bth 
partial  momenta  of  the  number  of  particles  in  the  size  class;  and  K  is  the  instmment  or  system  matrix  whereby 
element  represents  the  diffraction  contribution  of  a  unit  measure  of  particles  in  the  size  class  onto  the  ith 
detector.  Tne  elements  of  K  are  given  by: 
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where  WQ((x,aj)  is  a  weighting  function  for  the  jth  size  class  and  ky(a,0p  is  a  getieral  scattering  function  which 
gives  the  scattering  contribution  of  a  unit  quantity  of  particles  of  size  a  into  angle  0.  The  solution  or  measured 
particle  size  distribution  indicated  by  N  in  Eq.  (2)  can  in  theory  be  obtained  by  inverting  the  matrix  K. 

EXPERIMENTAL  PR0C:EDURE 

Hg.  1  shows  a  schematic  for  a  laser  diffraction  particle  sizing  instrumenL  A  laser  beam  is  spatially  filtered, 
eiqianded  aixl  collimated  to  a  few  millimeter  1/e^  diameter.  Particles  within  the  collimated  bem  scatter  light  which 
is  collected  by  a  receiving  lens  and  refracted  onto  the  detector  jriane.  A  commercially  availaUe  Insitec  particle 
sizing  detector  bead  and  electronics  were  used  for  the  eiqieriment .  The  layout  of  the  Insitec  detector  is  well-suited 
for  examining  detector  response,  given  its  odd  numbered  rings  on  one  half  of  the  detector  and  even  numbered  rings 
on  the  other.  This  design  has  large  non-iesponsive  sectiras  between  adjacent  rings,  which  provides  good  isolation 
between  adjacent  ctetectots. 

The  Insitec  detector  was  mounted  in  a  cantilevered  position  on  a  programmable  x-y  translation  stage  which 
moved  in  a  friane  normal  to  the  beam  axis.  A  HeNe  laser  beam  was  focussed  to  a  Utr  spot  25pm  in  diameter  and 
the  spm  was  positioned  in  the  small  bole  in  the  center  of  the  detector.  (The  bole  allows  the  strong  unscattered  light 
at  the  center  of  the  detector  to  pass  through,  thereby  reducing  stray  reflections  and  inner  ring  cross-taOc).  Coarse 
centering  was  performed  by  moving  the  detector  until  the  light  reflecting  oH  the  region  surrounding  the  bole 
disqipeaied  (ideating  the  light  was  passing  through  the  bole).  The  beam  power  was  adjusted  so  that  when  the 
laser  was  positioned  in  the  middle  of  an  active  detector  the  resultant  signal  was  just  below  the  saturation  level 
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Fig.  t.  Schematic  of  ceavcnticnal  Imer  diflraetioa  particle  aizing 
ayatam. 


Pig.  2.  Groaa  acction  of  an  Inailec  photodiode  detector  ahowiiy  the 
ideal  raapeoaiva  region  (R-octive),  the  transition  regioo  (grp).  Mid 
the  non-active  icgioo  (meeked). 


Several  runs  were  made  by  scanning  across  the  entire  detector  on  a  line  passing  through  the  center  in  one  direction 
and  then  rotating  the  detector  90  degrees  aixl  repeating  the  scan  in  the  other  direction.  Signals  were  recorded  at 
1pm  increments  and  the  entire  eiqietimental  apparatus  was  shielded  from  stray  room  li^L 

Before  scanning  the  detector,  the  rings  were  examined  under  a  microscqie.  Three  distinct  regions  were  noticed 
which  are  shown  in  Hg.  2  and  are  labeled  as  follows:  the  ideal  reqmnsive  or  active  region  (Rgctive)*  ^  transitioa 
re^on  (gap),  and  the  non-acdve  region  (masked).  Bounding  each  detector  is  a  thin  aluminum  film,  whidi  has  been 
deposit^  to  cover  most  of  the  area  between  detectors.  It  is  assumed  that  the  metal  helps  block  radiation  so  that 
masked  regions  are  not  reqioosive.  Though  Rg.  2  shows  just  one  detector,  all  detectors  were  examined  and  found 
to  be  similar. 

To  measure  the  gap  widdi,  the  detector  was  placed  on  a  two-axis  translation  stage  that  was  coupled  to  LED 
disfdays  showing  relative  stage  position.  The  resolution  of  the  display  was  1pm.  A  camera  connected  to  a  19-indi 
TV  screen  was  mounted  on  the  microscope  wUch  allowed  for  easy  visual  inspection  of  the  object  A  small  piece  of 
Ufte  was  placed  on  the  screen  to  mark  the  starting  and  ending  points  for  various  features.  An  inner  gap  edge  was 
visually  lined  up  with  an  imaginaiy  tangent  line  on  the  edge  that  would  be  parallel  to  the  tape.  The  stages  were 
moved  until  the  outer  gip  edge  was  reached  and  the  difference  between  stage  positions  was  recorded.  In  this 
manner  the  width  of  each  gap  (both  inner  and  outer  per  detector),  was  determined. 


RESULTS  AND  DISCUSSION 


The  average  gap  width  was  approximately  8pm  (results  from  all  detector  rings  can  be  found  in  Table  1).  The 
largest  measured  gq>  width  was  10mm  while  the  smallest  gt^  width  was  6mm,  though  most  widths  were  Smm.  The 
table  begins  with  detector  5  because  the  first  four  inner  tings  were  damaged  when  a  hole  was  laser  drilled  through 
the  detector  center.  Also  included  in  Table  1  is  die  visual  measured  width  of  the  Insitec  detectois  (from 
microscope). 

Rg.  3  displ^rs  the  individual  ring  signal  vs.  radial  distance  for  the  first  five  inner  odd-numbered  rings.  The 
rectangular  boxM  represent  the  specified  width  of  the  detectors  and  therefore  die  ideal  top  bat  leqionse  distributitm, 
while  the  curves  lefnesent  the  measured  signal.  The  measured  signal  is  a  convolution  of  the  detector  qiatial 
response  and  the  Gaussian  intensity  profile  of  the  laser  beam.  Since  the  measured  curves  are  centered  around  the 
sp^fied  poitHs,  it  was  concluded  timt  the  laser  beam  was  qiproximately  centered  ihout  the  detector.  The 
maximum  sigrud  is  dose  to  8  vdts,  however  diis  value  is  not  quite  reached  on  rings  4  and  S  because  the  width  of  the 
rings  is  less  than  twice  the  beam  width,  (where  twice  the  beam  width  represents  more  than  99%  of  the  total  beam 
energy).  Also,  overl^i  between  rings  occur  because  of  the  beam  width  and  ed^  effects  resulting  from  electraa 
bole  pair  and  photon  diffusion.  Note  bow  overiap  decreases  with  increasing  ring  size,  ududi  is  due  to  the  larger 
regions  of  non-reqionstve  area  between  rings.  Tbe  curves  appear  smooth  excqit  for  the  saturation  region  where 
small  deviations  in  signal  are  noticed.  It  should  be  mention^  that  KXK)  signals  are  avenged  at  eadipositioa 

In  order  to  be  able  to  calculate  wg  for  use  in  Eq.  (1)  it  is  necessary  to  model  tiw  gap  reqiooae.  This  reqionse  can 
be  estimated  by  deconvolving  the  meamred  dgnal  the  known  Oausdan  imendty  profile  of  the  laser  beam. 
After  several  runs  were  taken  of  the  inner  odd  rinp  (results  of  tings  11,13,15)  an  avenge  teqronse  in  the  gap  tegioo 
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Pig.  3.  Plot  efdeieetor  signal  VB.radia]  distance  for  five  famer  odd 
rings.  The  raelanglbe  rapresenl  the  apecified  dateeior  diroenaioiia. 
The  carves  are  the  aignak  reaalling  frein  a  \/tr  2StUn  diameter 
laser  beam  ct  wavolmgih  0.6328)iin.  Measaremenla  were  made  at 
Iftm  ineramenls. 


to  form  one  carve. 


average 
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wag  calculated  by  detennining  the  avenge  signal  for  each  ring  id  each  Ifun  position  (over  a  l(X)pin  distance)  and 
then  superimposing  the  data  to  calculate  an  overall  avenge.  This  average  gap  leqwnse  curve  is  shown  in  Hg.  4. 

Initial  insight  into  a  possible  model  for  die  gap  response  is  found  by  looking  for  symmetry  within  Hg.  4.  The 
symmetry  present  in  the  curve  points  to  a  potential  model  based  on  a  simple  fonedon,  most  likely  e:q>onemial.  A 
logical  model  to  predict  the  gi^  reqionse  based  on  physical  conditions  would  include  uniform  re^nse  over  the 
idealized  d^ector  region,  decaying  eiqionential  response  over  the  8pm  gap  width  and  no  reqionse  in  the  masked  or 
noo>ieq)onsive  regioa 

Photons  incident  mi  die  gap  regioa  art  absorbed  by  the  silicon  vriiere  they  create  eiectron-hede  pairs.  The 
electron-hole  pairs  begin  difliising  through  die  material  until  they  either  readi  the  p-n  juncdmi  wbm  there  is  a 
possibility  that  they  will  contribute  to  the  overall  cutim,  or  they  lose  their  energy  thrwgb  collisions  and  reemnbine. 
For  an  overview  of  the  physics  of  photodiode  operation  see  Sze  [8].  As  far  as  a  model  is  concerned  only  one 
parameter  will  need  to  determined;  the  "diffusion  length  constant"  governing  the  decaying  eiqionent^ 

However,  it  should  be  mentioned  that  when  convolving  the  Gaussian  beam  with  the  giqi  response  function  the 
value  of  the  dififosion  length  constant  has  only  a  slight  impact  on  the  slope  of  the  curve,  thereby  making  it  difiScult 
to  esdmate  the  proper  value  of  the  constant  to  match  the  measured  data.  Rg.  5  confirms  this  t^servadon,  where 
three  carves  rqiteseitting  the  same  beam  diameter,  but  different  difiiision  length  constants  are  displayed.  As  can  be 
imagined  the  gip  reiqionse  resulting  from  different  diffusion  constants  is  highly  variable,  yet  because  the  beam 
width  is  much  greater  than  the  gap  width  the  slopes  of  the  corves  are  nearly  equal.  If  the  beam  width  was  less  than 
the  gap  width  then  the  dopes  of  the  three  corves  would  be  much  different  Rirthermoie,  if  the  beam  was  reduced  to 
a  point  source  then  the  reqionse  would  simply  be  the  exponentiai  itself.  However,  even  with  large  beam  widths  a 


S.  Plot  of  dieoratical  carve*  calealatad  by  canvolvin|  a 
aa  beam  diameter  arilh  me  dieeialieal  detector  leapotiae.  The 
apectfied  detector  width  ie  SSfUn  which  ceneeponda  to  riiM  9  of 
the  iaailee  detector.  Calealation*  are  made  adlh  different  diffiuioo 
length  eontlaala  aa  indicated.  The  ideal  active,  gep,  and  maaked 
regkn*  are  alao  indicated  on  the  figate. 


Pig.  6.  Compariaen  of  theoretical  and  meaaared  aignal  va.  relative 

diatance  for  ring  9.  The  theotetieal  aignal  waa  ctJcalaled  baaad  on 
a  2jpm  beam  diameter,  l|tm  edge  width,  90|im  dffaaien  length, 
and  SSflm  detector  width. 


diffeienoe  exists  between  the  curves.  The  results  are  sensitive  to  difliision  length  in  that  the  curves  eiqtand  with 
inoeasiag  diCfiision  length  constant. 

Since  the  specified  width  of  each  detector  is  known,  a  theoretical  calculation  could  use  diffusion  constants  as  a 
means  to  expand  or  contract  the  curve  to  match  the  measured  data.  First,  the  total  laser  power  in  the  Gaussian  beam 
is  normalized  to  one  and  the  beam  is  moved  in  one  micron  increments  over  die  idealized,  gq>,  and  madced  regions. 
At  each  position  the  eiqiecied  signal  is  calculated  by  convolving  the  response  function  (assuming  a  diffusion  length) 
with  the  Gaussian  beam  intensity.  (As  can  be  seen  from  Rg.  2  there  are  instances  when  the  width  of  the  beam 
coven  all  three  regions  of  the  detector.)  These  values  are  then  compared  to  the  measured  dau  to  find  the  diffusion 
length  constam  that  provides  the  best  fit 

Eiqierimental  dau  from  ting  9  were  compared  to  theoretical  data  calculated  based  on  different  diffusion  length 
constanu.  Assumed  values  in  the  theoreticid  calculations  include  a  25itm  1/e^  Gaussian  intensity  profile  diameter 
laser  beam,  gap  width  of  8pm,  and  total  deteaor  width  obtained  from  visual  inspection.  Results  stowed  that  a 
diffiisioo  length  constant  of  50pm±10pm  (which  falls  within  observed  measurements  [8])  gave  the  best  fit  between 
experimental  and  theoretical  values.  A  plot  of  ting  9  normalized  sigiul  vs.  distance  is  found  in  Hg.  6,  and  as  can  be 
seen  from  the  grifih  the  fit  is  quite  good.  As  a  test  for  the  motol,  comparisons  were  made  between  erqierimental 
and  measured  resets  for  many  of  the  other  tings  and  similar  results  were  noted.  Thus,  it  was  conchitM  that  a 
simple  exponential  function  ^ipeats  to  adequately  predict  the  giq>  resfioase  of  the  Insitec  ring  detector.  The 
weighting  function  over  the  gap  width  can  to  written  as: 

^dedge  *  p/p  “  W 

where  p(6)  is  the  local  responsivity,  p  is  the  maximum  responsivity  in  the  ideal  active  region  of  the  detector,  xd  is 
the  difiusion  length  constant,  and  x  is  the  position  measured  from  the  ideal  edge  into  the  gq>. 

Now  that  the  weighting  function  wq  has  been  determined  a  new  instrument  matrix  K  can  to  calculated  by 
substituting  the  measured  gap  weighting  function  of  Eq.  (4)  into  Eq.  (3)  in  place  of  tiie  usually  assumed  step 
functiotL  To  quantify  the  eSeet  of  the  detector  wei^ting  function  on  the  instrumem  matrix,  two  matrices  were 
calculated;  one  assuming  a  step  detector  reqionse,  the  other  the  actual  measured  giq>  response  (and  weighting  values 
of  0  and  1  in  the  masked  and  idealized  regions,  respectively).  A  plot  of  the  difference  matrix  (elemetk  by  element) 
between  the  two  matrices  is  presented  in  Hg.  7.  Tto  matrices  were  calculated  on  an  area  basis  with  rectsogular 
weighting  functions,  uniform  1^  volume  withm-ciass  distributions,  and  a  Fraunhofer  diffraction  approximation. 

The  surface  plot  shows  almost  no  change  on  the  outer  rings,  while  the  difference  increases  as  ring  geometry 
decreases.  This  is  eiqiecied,  as  tto  ratio  of  giq>  area  to  visual  detector  area  is  much  greater  for  smaller  rings  than 
larger  ones. 
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Pig.  7.  Plot  of  the  difloTBOM  matrix  between  K  baaed  on  a  alep 
leaponae  detector  weighting  Ametion  and  K  obtaiiied  with 
dworetieal  detector  leaponae  weighting  fcnetion  (which  eecoonte 
for  detector  reeponae  fai  the  gap  legion  between  Ite  meak  and 
Ueahsad  detector  boundary). 


Pig.t.  Plot  of  calibratioa  faelore  va.  detector  width  for  varioai 
inainNiienU  and  calibinlioa  lechnitpiet.  The  numbeca  in  the  legend 
corieepond  to  the  following:  I)  Malvem  2d00,  inetrwincnt  #1, 
uniform  ineoheroil  illumination  (6,1 1).  2)  Malvem  2600, 
inilrBmentff2,factotyBoppt>ed^!bfelionfaemra(91.  3)  Malvem 
2600,  iaatrument  #3,  uniform  bieoheient  illumination  [10].  4) 
Inahee,  unHoim  incoheieat  iHumination.  3)  Inailee,  Iheoreticd 
piedictient  beaed  on  gap  leaponte  method. 


It  bag  already  been  meatiooed  that  it  is  possible  to  predict  calibiatioo  factors  based  on  gq>  or  traasition  region 
response  as  opposed  to  systematic  betweeo-detector  vaiiadoos  in  ibe  lesponsivity.  Tbe  accuracy  of  tbe  "g^" 
calculated  faaon  is  based  on  tbe  piemiae  that  tbe  major  contribution  to  calibration  factors  is  from  gap  response.  A 
predicted  calibration  factor  is  defined  bere  as  tbe  rMio  of  tbe  ideal  detector  ou^t  to  tbe  actual  detectm  output  when 
illuminated  widi  light  of  uniform  intensity.  The  ideal  detector  output  is  calculated  by  iraegraling  tbe  iix:id«t  light 
intensity  over  the  ideal,  visual  detector  area  vdiile  assuming  that  tbe  weighting  function  W9  is  constant  and  equal  to 
one.  Tbe  actual  detector  ouqiut  is  found  in  tbe  same  manner  except  tbe  area  over 'vhich  the  integration  is  petfonned 
includes  tbe  giqi  region  and  tbe  associated  wei^itiog  function  In  equation  form  we  have: 


^pred 


/idA 


(5) 


where  Cpj^  is  tbe  predicted  calibration  factor  and  I  is  die  intensity  (W/m^. 

Tradinoiially,  calibration  factors  have  been  determined  based  on  illuminating  tbe  detectors  with  uniform  li^t 
intenstty.  We  determined  two  sets  of  cafibration  factors  for  tbe  Insrtec  detector.  One  i  of  factors  was  calculated 
using  Eq.  (5),  while  tbe  other  set  was  determined  using  tbe  uniform  light  metbod.  Rr . ^Its  of  tbe  two  methods  are 
disfdayed  in  Rg.  10  along  with  calibration  factors  from  3  different  [9>11]  Malvem  2600  instruments  determined 
empirically  with  tbe  uniform  light  method.  Since  tbe  Insitec  detector  geometry  differs  sli^tly  from  the  Malvem 
detector  geometry,  tbe  calibration  frKtots  are  plotted  vs.  detector  width.  As  is  seen  from  tbe  plot  the  curves  are  in 
reasonable  agreemenL  Noticed  trends  are  that  all  factors  decrease  in  value  as  detector  area  decreases. 

Consequently,  tbe  inner  detectors  over*tespond  relative  to  tbe  outer  dmectois.  Calibratioa  factors  determined  from 
Eq.  (5)  are  limited  by  tbe  fact  that  experimental  data  wme  collected  along  two  lines  tunning  through  tbe  center  of 
tbe  detector  instead  of  scamting  tbe  entire  detector  surface.  Additionally,  in  Eq.  (5)  it  was  assumed  that  tbe 
weigtuing  functioo  equals  a  constam  over  tbe  ideal  detector  region,  however,  Hg.  4  shows  slight  variation  in  tiiis 
regioa  One  other  area  ofuncertainty  which  we  did  not  investigate  is  tbe  ends  of  each  detector.  Though  we 
accounted  for  tbe  aid  regioa  in  tbe  calculations  tbe  regioa  a  (tifferent  from  tbe  rest  of  tbe  perimeter  regions  in  that  a 
masked  mea  is  iMt  found  8|tm  from  tbe  ideal  detector  boundary.  A  drawback  to  the  uniform  light  metbod  is  that 
linear  detector  response  is  assumed  with  no  cross-talk.  If  either  of  tbe  above  assumptions  is  not  correct  then  errors 
in  calibcttioa  factors  will  propagate  through  tbe  remaming  factors  that  need  to  be  determined.  This  results  from  tbe 
requirement  tint  values  acquired  from  differem  light  levels  must  be  qUioed  togetiier  Pj.  Overall,  it  ^ipears  that 
camnation  fiictofs  calculate  using  Eq.  (5)  compare  reasonably  well  with  factors  determined  with  tbe  uniform  li^ 
m^wd.  especially  conndering  tbe  fact  that  factors  are  fm  different  daecters  from  different  laboratories. 
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CONCLUSIONS 


In  conclusion,  the  transition  region  response  has  been  measured  and  modeled  for  an  Insitec  detector.  A  simple 
exponential  with  a  diffusion  length  constant  of  SOpm  accurately  predicted  the  experimental  Haf  The  detector 
weighting  function  has  been  calculated  and  incorporated  into  the  instrument  matrix  K.  Calibration  factors  for  the 
Insitec  detector  have  been  predicted  based  on  edge  effects  and  were  compared  with  calibration  factors  derived  from 
the  uniform  light  source  technique.  Both  methods  indicated  that  the  inner  rings  over-respond  relative  to  the  outer, 
larger  rings.  It  is  difScult  to  conclude  which  method  produces  the  mote  accurate  results.  The  uniform  light  method 
has  the  disadvantage  that  detector  signals  resulting  from  several  light  levels  must  be  measured  atxl  the  calibration 
curves  spliced  together.  The  gap  response  model  has  the  disadvaruage  that  it  is  based  on  2-dimensional  results  and 
it  must  be  assumed  that  gap  response  is  axisymmetric.  Another  independent  technique  for  mlmlnting  calibration 
factors  could  help  verify  whether  the  gap  or  uniform  light  source  is  more  accurate.  The  effect  of  using  the 
theoretical  calibration  factors  instead  of  the  uniform  light  source  factors  would,  in  general,  translate  into  slightly 
larger  particle  sizes  because  of  the  higher  signals  on  the  inner  rings. 
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General  solution  to  the  inverse 
near-forward-scattering  particle-sizing  problem  in 
multiple-scattering  environments:  theory 


E.  Dan  Hirleman 


A  general  solution  to  the  problem  of  measuring  the  size  distribution  of  large  particles  in  optically  thick 
media  by  using  small-angle  li^t  scattering  is  presented.  The  qjproach  is  general  in  the  sense  that  no 
assumption  of  the  form  of  the  particle-size  distribution  function  is  required,  although  the  particles  must 
be  distributed  uniformly  throu^out  the  medium.  The  method  is  based  on  a  successive  order,  discrete 
ordinates  approach  for  modeling  multiple-scattering  phenomena  and  requires  that  the  particle  field  be 
interrogated  by  using  an  array  of  near-forward  input  light  angles.  The  scattering  redistribution  matrix  is 
thereby  determined,  which  permits  a  numerical  inversion  of  the  problem  to  obtain  the  single-scattering 
signature.  Finally,  conventional  inverse  scattering  methods  are  used  to  reconstruct  the  partkle-size 
distribution  from  the  near-forward  (sin^e-scattering)  light-scattering  pattern. 


Introduction 

The  phenomenon  of  multiple  scattering  significantly 
complicates  the  analysis  of  radiation  transfer  throu|^ 
particulate-laden  m^a.  Nonintrusive  diagnostics  for 
various  flows  depend  on  an  understanding  of  light 
scattering  and  propagation  and  for  this  reason  can 
unfortunately  not  be  used  in  the  many  important 
systems  and  applications  that  involve  optically  thick 
media.  We  develop  a  general  solution  to  the  problem 
of  measuring  the  size  distribution  of  large  particles  in 
opticeiUy  thick  media  using  small-angle  li^t  scatter¬ 
ing.  Numerical  experiments  are  used  to  demonstrate 
the  feasibility  of  the  technique. 

In  a  previous  paper  by  Hirleman'  a  theoretical 
approach  to  solving  the  direct  problem  (i.e.,  calcula¬ 
tion  of  the  scattering  signature  given  the  particle  size 
distribution)  for  the  case  of  near-forward  scattering 
by  particles  that  are  large  compared  to  the  wave- 
len^h  in  optically  thick  media  was  presented.  While 
other  models  exist  for  predicting  the  scattering  prop¬ 
erties  of  thick  media,  the  discrete  ordinates,  succes¬ 
sive  order  approach  of  this  previous  work'  was  devel¬ 
oped  specifically  for  incorporation  into  schemes  for 
solving  the  inverse  problem.  We  present  a  general 
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solution  for  the  corresponding  inverse  problem,  i.e., 
that  of  estimating  the  size  distribution  function  of  an 
ensemble  of  large  particles  by  using  mathematical 
inversion  of  the  near-forward-optical-scattering  prop¬ 
erties  of  the  optically  thick  medium. 

Theory 

Consider  a  particle-laden  medium  illuminated  by  a 
planar  electromagnetic  wave  at  an  arbitrary  (but 
small)  angle  relative  to  the  laboratory  optical  Z 
axis.  The  coordinate  system  of  interest  is  shown  in 
Fig.  1  with  a  sin^e  p^icle  at  the  origin  and  for  a 
wave  incident  at  =  0.  Assume  that  the  scattering 
phenomenon  is  axisymmetric,  i.e.,  with  no  azimuthal 
<I>  dependence  of  the  scattered  field  as  would  be  the 
case  for  an  effectively  infinite  medium  of  spherical 
particles.  If  the  incident  field  actually  consist^  of  the 
superposition  of  waves  covering  all  possible  values  of 
Ojnj  for  a  given  (i.e.,  a  hollow  cone  of  incident 
li^t),  the  overall  scattering  characteristics  of  the 
medium  would  be  ssmunetric  about  the  optical  axis.  If 
nonspherical  particles  were  present,  so  that  the  indi¬ 
vidual  scattering  signatures  were  not  axisynunetric, 
the  overall  radiation  transfer  problem  could  still  be 
considered  axisymmetric  if  there  were  a  large  number 
of  randomly  oriented  particles  in  the  me^um.  It  is 
also  assfuned  that  the  scattering  processes  involved 
are  linear  (in  the  sense  of  a  linear  S3rgtem).  The  linear 
system  assumption  require  that,  if  the  incident 
optical  power  is  multipfied  by  a  constant,  the  scat¬ 
tered  energy  leaving  the  medium  in  all  dii^ions  is 
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Fig.  1 .  System  for  analysis  of  multiple  scattering.  The  incident  electromagnetic  energy  is  shown  as  a  photon  passing  through  the  origin  of 
the  laboratory  XY2  coordinate  tystem  traveling  in  a  direction  defined  by  angles  and  A  plane  wave  traveling  in  the  same  direction 
would  be  scattered  by  the  spherical  particle  shown  at  the  origin  in  a  scattering  pattern  axisymmetric  with  respect  to  the  Poynting  vector, 
i.e.,  axitymmetric  about  the  z'  axis.  The  scattered  light,  indicated  as  a  scattered  photon  in  the  figure,  is  transformed  back  into  the  laboratory 
coordinate  system  before  analysis  of  the  next  scattering  event. 


similarly  multiplied  by  the  same  factor.  Finally,  depen¬ 
dent  scattering  is  neglected;  i.e.,  it  is  assvimed  that 
scattered  fields  add  incoherently  (on  an  intensity 
rather  than  amplitude  basis).  It  is  important  to  note 
here  that  the  total  scattered  energy  that  leaves  the 
medium  in  any  direction  is  not  necesborily  propor¬ 
tional  to  the  number  of  particles  in  the  system;  i.e., 
the  energy  scattered  in  a  given  direction  by  an 
Hf-particle  medium  is  not  equal  to  M  times  the 
scattering  contribution  of  one  isolated  particle  in 
conditions  of  significant  multiple  scattering.  How¬ 
ever,  the  classic^  definition  of  a  linear  system  in¬ 
volves  proportionality  between  input  and  output  vari¬ 
ables  (in  this  case  the  incident  and  scattered  energy) 
but  does  not  constrain  the  effects  of  changing,  say, 
the  extent  or  population  density  of  the  medium. 

The  approach  used  here  is  based  on  d  ’^tizingthe 
range  of  incident  and  scattering  angk  nto  a  finite 
number  of  dir  Ttions,  each  representing  a  finite  range 
of  scattering  angles.  This  term  discrete  ordinates  has 
been  used  by  van  de  Hulst^  and  others  to  describe  this 
general  method.  For  the  axis}rmmetric  problem  of 
interest  here  the  discrete  scattering-an^e  ranges  are 
independent  of  the  azimuthal  scattering  an^e  <I>  and 
subtend  2ir  rad  in  for  discrete  values  of  0.  These 
discrete  values  of  incident  and  scattering  cones  repre¬ 
sented  1^  6  are  then  the  discrete  or^nates  of  the 
system.  We  define^  a  scattering  vector  S  so  that  the 
ith  element  of  S  contains  a  measure  of  the  radiant 


energy  proceeding  in  directions  represented  by  the 
angle  6^.  The  effect  of  a  scattering  medium  then  is  to 
operate  on  the  incident  light  and  redistribute  the 
radiant  energy  over  the  various  0^.  (Some  of  the  light 
may  also  be  absorbed  in  the  medium,  and  this  effect  is 
considered  below.)  It  is  convenient  to  define  a  scatter¬ 
ing  redistribution  matrix  H„,  so  that  the  matrix 
element  in  the  ith  row  amd  the yth  column  H,„(i,  j)  is 
the  gain  or  efficiency  with  which  optical  energy 
incident  in  a  direction  represented  by  0^  is  redistrib¬ 
uted  by  the  medium  into  direction  0,  .  The  effect  of  the 
linear  medium  can  be  written  as 

(1) 

where  Si„.  is  a  vector  so  that  theyth  element  repre¬ 
sents  the  optical  energy  incident  on  the  medium  in 
direction  0,  and  S„  is  a  vector  so  that  the  ith  element 
represents  the  amount  of  optical  energy  exiting  the 
medium  in  direction  0,.  Note  that  Sj„c  ^ 

somewhat  analogous  to  the  Stokes  vectors  from  the 
electromagnetic  scattering  theory,*  except  that  polar¬ 
ization  is  not  considered.  The  assiunption  that  polar¬ 
ization  effects  are  ne^igible  is  reasonable  for  the 
near-forward-scattering  regime  of  interest  here.  The 
subscripts  m  on  and  S„  indicate  that  these  iq>p]y 
to  the  general  multiple-scattering  situation,  i.e.,  1^ 
represents  the  redistribution  effect  of  a  general  me¬ 
dium  and  S„  is  the  scattering  signature  that  may 
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include  multiple-scattering  effects.  These  are  distinct 
from  a  single-scattering  redistribution  matrix  H,  and 
scattering-order  vectors  S„  are  defined  below.  is  a 
function  of  the  particle-size  distribution  and  the 
optical  depth  of  the  medium. 

Selection  of  the  appropriate  number  of  and  the 
specific  values  of  the  ^screte  scattering  angles  0^  (i.e., 
the  discrete  ordinates)  depends  on  the  scattering 
properties  of  the  medium.  The  number  of  angles 
considered  is  an  important  factor  that  controls  the 
accuracy  with  which  the  discrete  approximation  in 
Eq.  (1)  approximates  the  actual  radiation  transfer 
process  that  is  actually  continuous  in  0  and  d>.  We 
select  0i  to  be  consistent  with  the  photodetector 
geometries  used  in  our  experimental  work  on  particle 
sizing.  It  has  been  shown  by  Hirleman®  that  the 
optimal  detector  geometry  for  particle  sizing  in  the 
large-particle,  small-scattering  angle  regime  (where 
Fraunhofer  diffraction  theoiy  is  valid)  consists  of 
log-scaled  rings.  In  the  theoretical  work  reported  here 
we  use  collection  aperture  geometries  corresponding 
to  the  ring  elements  of  the  RSI  ring/wedge  detector'* 
for  which  dimensional  data  are  available.^  We  con¬ 
sider  optical  energy  proceeding  nominally  along  the 
optical  axis  as  the  ^st  element  of  and  S„„  so  that 
this  element,  defined  by  the  representative  angle  0o, 
subtends  a  range  of  scattering  angles  from  zero  to  a 
small  angle  corresponding  to  the  central,  circular 
detector  element  used  for  measuring  extinction.  The 
angle  range  considered,  i.e.,  @o  -♦  0„,  where  n  +  1  is 
the  number  of  discrete  ordinates,  must  be  sxifficiently 
large  that  only  a  negligible  amount  of  optical  energy 
falls  outside  it.  For  that  reason  the  value  of  0„ 
required  for  a  particular  problem  depends  on  the 
scattering  properties  of  the  medium  (i.e.,  through  the 
peuticle-size  distribution)  and  the  optical  depth  of  the 
medium. 

The  conventional  method  for  particle  sizing  by 
using  near-forward  Ught-scattering  signatures  in¬ 
volves*®®  interrogation  of  the  medium  with  a  colli¬ 
mated  input  beam  that  defines  an  optical  axis  (0  =  0) 
that  is  centered  in  the  ring  detector.  In  this  case  the 
normahzed  incident  energy  vector  8;,,^  is  given  by 

S,,,  =  [1,  0,  0, . . .  of.  (2) 

Now  from  Eq.  (1)  we  see  that  in  the  conventional 
system,  with  the  angular  distribution  of  the  incident 
optical  energy  given  by  Eq.  (2),  the  output  scattering 
signature  S„  is  just  equal  to  the  first  column  of  H„,. 
The  direct  problem,  i.e.,  calculation  of  the  output 
scattering  signature  given  an  input  beam  and  a 
medium  defined  by  the  particle-size  distribution  tmd 
an  optical  depth  b,  can  Ite  viewed  as  determining  H,,, 
and  using  Eq.  (1)  to  predict  S„.  The  inverse¬ 
scattering  particle-sizing  problem  is  then  a  matter  of 
determining  a  particle-size  distribution,  which  in 
theory  should  have  produced  the  observed  scattering 
signature  S„.  To  do  this  it  is  in  turn  necessary  to  be 
able  to  calculate  theoretically  the  matrix  H,,,  vdid  for 


an  arbitrary  particle-size  distribution  and  optical 
depth. 

If  it  is  desired  to  extract  n  unique  pieces  of  informa¬ 
tion  about  the  particle-size  distribution  in  the  me¬ 
dium,  it  is  necessary  to  make  at  least  n  measure¬ 
ments.  Often  the  ra  pieces  of  information  (unknowns) 
are  the  quantity  of  particles  in  n  discrete  size  classes, 
and  the  measurements  of  the  h^t  scattered  at  n 
discrete  angles  provides  n  equations  to  produce  an 
ra  X  n  Unear  system.  Now  interrogation  of  a  medium 
with  an  incident  field  as  indicated  by  Eq.  (2)  resulted 
in  an  experiment  that  produced  n  measurements,  i.e., 
the  elements  of  S„.  Further  consideration  reveals 
that  Eq.  (1)  can  in  fact  support  ra  independent  experi¬ 
ments,  each  consisting  of  scattering  measurements  at 
ra  angles  (i.e.,  ra®  total  independent  measurements). 
Any  of  these  experiments  could  be  used  to  obtain  an 
estimate  of  the  particle-size  distribution  function, 
although  conventional  near-forward-scattering  meth¬ 
ods  use  only  one.  An  vector  with  the  yth  element 
as  the  only  nonzero  element  represents  a  hoUow  cone 
of  illumination  with  half-angle  0^  .  For  an  experiment 
of  this  type  only  the  j  th  column  of  H„  would  enter  the 
problem.  An  optical  system  designed  to  produce  hol¬ 
low-cone  illumination  as  discussed  is  shown  in  Fig.  2. 
A  programmable  mask  with  an  array  of  ra  annular 
shutters  is  placed  in  the  front  focal  plane  of  a 
transmitting  lens.  This  ring  or  cylinder  of  Ught 
transmitted  through  the  mask  is  transformed  by  the 
lens  to  contain  only  a  small  range  of  angles  in  the 
back  focal  plane  of  the  lens;  i.e.,  a  hoUow  cone  with 
half-angle  0,  where  0  is  related  to  the  mean  radius  r 
of  the  annular  shutter  and  the  lens  focal  length  fhy 
0  =  rif.  A  matched  annular  shutter  array  is  used  at 


Fig.  2.  Simplified  schematic  of  an  optical  system  that  allows 
multian^e  interrogation  and  multian^e  scattering  measurements 
of  a  spray  or  particle  field.  A  programmable  annular  shutter  or 
mask  in  the  front  focal  plane  of  the  transmission  lens  creates  an 
annular  ring  or  (^linder  of  light.  This  is  then  converted  into  a 
hollow  cone  of  light  (all  light  traveling  with  nominally  the  same 
an^e  relative  to  the  optical  axis)  by  the  Fourier  transform 
process  of  the  lens.  Ldg^t  scattered  the  medium  is  transformed 
from  an  angular  to  a  spatial  pattern  at  the  back  focal  plane  of  the 
receiving  lens.  Li^t  scattered  into  discrete  ranges  of  <|>  and  6  are 
sampled  by  another  annular  shutter  array  at  the  detector  plane 
and  transmitted  to  a  detector  (not  shown)  behind  the  shutter 
array. 
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the  back  focal  plane  of  the  receiving  lens,  and  for  each 
incident  cone  angle  the  system  is  cycled  through  all 
rings  at  the  detector  plane.^ 

Theory  for  Determination  of 

It  has  been  pointed  out  that  Eq.  (1)  can  theoretically 
be  used  to  determine  particle-size  distributions  in 
arbitrary  media,  if  we  assume  that  the  elements  of  H,„ 
can  be  theoretically  expressed  as  a  function  of  the 
quantity  of  particles  in  n  discrete  size  classes.  How¬ 
ever,  the  sttjbility  of  the  proposed  inverse  solution  is 
in  question,  and  methods  for  calculating  are  not  in 
the  hterature.  In  this  section  a  derivation  of  is 
presented. 

Recall  that  H,„  is  defined  so  that  the  element  {i,j)  is 
related  to  the  probabUity  that  Ught  incident  at  6^  will 
leave  the  medium  in  a  direction  specified  by  6^.  This 
probability  in  turn  depends  on  the  single-scattering 
phase  function  and  the  extent  of  the  medium.  If  the 
medium  is  optically  thin,  multiple  scattering  can  be 
ne^ected  and  the  problem  is  relatively  simple  as  the 
phase  function  can  be  used  directly.  The  successive 
orders  approach  used  here'  effectively  reduces  the 
problem  to  an  analysis  of  a  sequence  of  single¬ 
scattering  events,  each  of  which  is  easily  modeled.  It 
is  in  this  context  that  the  single-scattering  redistribu¬ 
tion  matrix  H  is  defined  so  that  represents  the 
probability  that  a  photon  incident  on  a  medium  at  an 
angle  representative  of  the  jth  detector  will  leave  the 
next  scattering  event  in  a  direction  corresponding  to 
the  ith  detector.  Note  that  H  depends  on  the  detector 
geometiy  and  the  scattering  characteristics  of  the 
medium,  which  is  a  function  of  the  unknown  particle- 
size  distribution  but  not  on  the  optical  depth  or 
extent  of  the  medium.  The  matrix  H  is  needed  for  use 
in  a  rectirsive  scattering-order  equation: 

S.  =  HS,.,,  (3) 

where  S„_i  emd  S,  are  the  scattering  signature  vectors 
of  light  (photons)  scattered  exactly  n  -  1  and  n  times, 
respectively.  A  special  case  of  Eq.  (3)  is  for  unscat¬ 
tered  (incident)  energy: 

S,  =  HS,,.  =  HS,.  (4) 

Inspection  of  Eq.  (4)  indicates  that  the  first  column 
of  H  is  just  the  normalized  scattering  signature, 
which  would  be  obtained  from  an  opticE^y  thin 
aerosol  illuminated  with  a  collimated  be^.  Now  H  is 
theoretically  determined  (for  a  given  size  distribu¬ 
tion)  by  using  an  analysis  that  combines  some  trans¬ 
formations  that  relate  light-scattering  and  laboratory 
coordinate  systems  with  the  calculation  of  scattering 
phase  functions.  The  development'  requires  a  small- 
angle  approximation  and  therefore  may  be  considered 
to  be  \^d  for  relatively  large  particles.  The  small- 
an^e  assumption  simplifies  ^culation  of  the  proba¬ 
bility  distribution  over  the  number  of  scattering 
events  that  are  related  to  the  physical  (as  opposed  to 
optical)  medium  width  throu^  the  optical  depth  b. 


Now  the  composite  scattering  signature  S„  (super¬ 
position  over  all  scattering  orders  So  „)  is  given  by' 

8.  =  exp(-6)[I  +  a^H  +  (a^f/2lf‘  +  +  ■  ■  -IS,,,  (5) 

where  I  is  the  identity  matrix  and  is  defined  as  a 
forward-scattering  albedo  where  =  0.5  in  the 
difiraction  regime  of  interest  here. 

Using  the  definition  for  the  logarithm  of  a  matrix  in 
combination  with  Eq.  (5)  we  obtain 

S„  =  exp(-6)exp(o^H)Si,,  (6) 

and  by  comparing  Eqs.  (1)  and  (6)  we  obtain 

H„  =  exp(-6)exp(a^H),  (7) 

which  provides  a  rather  straightforward  link  between 
the  single-scattering  redistribution  matrix  H  and  the 
general  medium  matrix  H„,.  Use  of  the  optical  system 
of  Fig.  2  permits  measurement  of  the  elements  of 
H,„  according  to  Eq.  (1).  Inverting  Eq.  (7)  we  obtain 

H  =  ln[exp(6)H„]/(Of6)  (8) 

which  provides  a  hnk  back  to  H  if  H,,,  cmd  b  are 
measui^  experimentally.  Recall  that  the  first  col¬ 
umn  of  H  is  just  the  single-scattering  signature, 
which  would  be  obtained  from  collimated  illumina¬ 
tion  of  an  optically  thin  medium.  Conventional  near¬ 
forward-scattering  instruments  routinely  use  this 
single-scattering  signature  for  0;,,^  =  0  in  an  inversion 
process  to  obtain  the  particle-size  distribution.  There¬ 
fore,  if  the  first  column  of  H  can  be  determined, 
conventional  single-scattering  inversion  schemes  can 
be  applied.  The  mathematical  procedure  given  in  Eq. 
(8)  whereby  operations  are  performed  on  an  experi¬ 
mentally  measured  H„,  is,  in  a  sense,  equivalent  to 
reaching  into  the  medium  to  collect  photons  that  have 
been  scattered  once  and  only  once. 

Relationship  to  the  Inverse  Single-Scattering  Problem 

The  single-scattering  inverse  problem  is  generaUy 
written®  in  a  form 

S  =  K  •  N,  (9) 

where  S  is  the  scattering  signature  over  the  detec¬ 
tors,  N  is  the  unknown  size  distribution  vector  so  that 
represents  the  quantity  of  particles  in  the  yth  size 
class,  and  K  is  the  scattering  matrix  so  that  element 
indicates  the  amount  of  energy  that  a  unit  quan¬ 
tity  of  particles  representative  of  the  yth  size  class 
would  scatter  into  the  ith  detector  assuming  single 
scattering.  Now  the  elements  of  K  can  be  determined 
by  using  the  Lorenz-Mie  scattering  theory  or  some 
approximation  such  as  the  Fraui^ofer  diffraction 
theoiy.  The  matrix  is  calculated  by  using  wei^ting 
functions  for  both  the  size  classes  and  the  detector 
geometries.® 

Now  it  would  be  useful  to  formulate  the  inverse 
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multiple-scattering  problem  as  a  lineeir  system  of  the 
form: 

S„  =  K„N,  (10) 

where  we  must  define  (and  determine)  as  the 
scattering  coefficient  matrix  that  applies  for  a  multi- 
ple-scatterii^  medium.  The  inverse  multiple-scatter¬ 
ing  problem  would  then  involve  inversion  of  the 
matrix  or,  in  the  case  of  an  ill-conditioned  K,., 
some  other  method  for  solving  the  linear  system 
given  by  Eq.  (10)  would  be  used.  Now  it  has  previ¬ 
ously  b^n  noted  that  the  matrix  H  depends  on  the 
particle-size  distribution  through  its  effect  on  the 
scattering  signature,  so  we  choose  to  define  a  dif¬ 
ferent  H  for  each  size  class  of  N.  If  we  denote  as  the 
single-scattering  redistribution  matrix  for  particles 
representative  of  the  yth  size  class  and  as  the 
multiple-scattering  redistribution  matrix  for  a  me¬ 
dium  containing  only  particles  in  size  classy,  Eq.  (1) 
becomes  for  this  case 

(11) 

where  is  the  multiple-scattering  signature  pro¬ 
duced  by  a  quantity  particles  from  size  classy  with 
illumination  specified  by  S,,„,  and  where  the  multipli¬ 
cation  by  the  scalar  Nj  is  a  scalar  product  but  the 
product  is  a  vector  operation.  When  we  ex¬ 
pand  Eq.  (11),  using  the  definition  in  Eq.  (7),  we 
obtain: 

=  exp(-<),)exp(yHo^^;StoJI;,  (12) 

where  yH  is  defined  as  the  single-scattering  redistribu¬ 
tion  matrix  in  a  medium  containing  only  particles 
representative  of  theyth  size  class  and  bj  in  Eq.  (12)  is 
defined  by 

b,  =  PjC^=(N,/V)C^/,  (13) 

where  p,  and  are  the  number  density  (number 
per  unit  volume)  and  mean  extinction  cross  section, 
respectively,  of  particles  in  size  classy , /is  the  physical 
length  of  the  medium,  and  V  is  the  volume  within 
which  the  number  quemtity  Nj  applies.  Note  that 
taking  the  particle  number  density  p,  =  (Nj/V)  re¬ 
quires  that  Nj  be  defined  on  a  number  basis  (as 
opposed  to  an  area  or  volume  basis).  The  quantities  p, 
and  would  be  obtained  by  integrating  over  the 
size  range  of  the y  th  size  class  by  using  an  appropriate 
within-class  wei^ting  function.^ 

To  return  to  a  form  such  as  in  Eq.  1  it  is  tempting  to 
use  a  superimposition  of  the  scattering  signatures  for 
the y  size  classes  from  Eq.  ( 12).  It  is  the  case  for  which 
the  composite-scattering  signature  could  be  calcu¬ 
lated  as  a  linear  combination  of  signatures  from  each 
class,  but  unfortunately  the  correct  scattering  signa¬ 
tures  to  be  used  in  the  summation  are  not  given  by 
Eq.  (12).  Equation  (12)  applies  to  a  medium  contain¬ 
ing  only  one  size  and  does  not  apply  in  the  general 
case  with  particles  from  a  mxiltiplicity  of  size  classes. 


This  nonlinearity  or  coupling  of  the  equations  would 
come  through  the  fact  that  the  redistribution  matrix 
for  a  given  size  class  in  a  polydisperse  medium 
depends  on  the  populations  of  all  the  other  size 
classes.  In  other  words  the  matrix  K,„  in  Eq.  (10) 
cannot  be  calculated  independent  of  a  knowledge  of 
the  particle-size  distribution  N,  and  hence  Eq.  (10)  is 
not  a  linear  system  in  the  conventional  sense.  How¬ 
ever,  the  matrix  K,„  could  be  calculated  for  a  given  N, 
and  the  predicted  scattering  signature  S  from  Eq. 
(10)  could  be  compared  with  a  measured  one  to 
provide  the  basis  for  an  iterative  solution  to  deter¬ 
mine  a  best  fit  N.  In  contrast  the  approach  suggested 
here  is  different.  The  rather  involved  calculations 
that  are  required  to  account  for  a  coupling  of  the 
multiple  scattering  between  size  classes  are  replaced 
here  by  a  number  of  independent  measurements  that 
serve  to  determine  this  coupling  experimentally.  The 
results  of  numerical  studies  of  the  performance  of  the 
new  approach  are  discussed  below. 

ResuHs  of  Numerical  Experiments 

To  determine  the  feasibility  of  the  proposed  approach 
numerical  experiments  were  performed  to  simulate 
the  process  of  (1)  measuring  b  and  H,„  by  using 
multiangle  interrogation,  (2)  calculating  an  estimate 
of  the  single-scattering  redistribution  matrix  H  by 
using  Eq.  (8),  and  finally  (3)  reconstructing  a  particle- 
size  distribution  function  by  inverting  Eq.  (9).  To 
calculate  H„,  it  is  first  necessary  to  obtain  a  theoreti¬ 
cal  H  to  use  in  Eq.  (7).  A  plot  of  the  calculated  matrix 
H  for  National  Bureau  of  Standards  SUM  1003a 
using  the  RSI  detector  geometry  within  an  optical 
system  with  a  receiving  (transform)  lens  focal  length 
of  125  mm  is  given  in  Fig.  3.  Note  that  the  axis  in  the 


(31.0) 


Fig.  3.  Plot  showing  the  calculated  sin^e-scattering  redistribu¬ 
tion  matrix  H.  The  redistribution  matrix  H  was  calculated  for  the 
particle-size  distribution  of  National  Bureau  of  Standards  (U.S.) 
SRM  1003A  for  a  laser  difiraction  system  that  uses  the  RSI 
detector  geometry*  ’  with  a  transform  lens  focal  length  of  125  mm 
amd  the  successive  order,  discrete  ordinates  multiple-scattering 
model.'  The  upper  left  element  of  H,  i.e.,  or  H(0,  0),  is  shown  on 
the  left  and  applies  to  li^t  incident  at  the  smallest  angle  6,  (0,  =  0 
or  collimated  input),  which  is  scattered  back  into  the  same  an^e 
(identically  forward  scattering)  by  the  next  scattering  event.  The 
lower  point  in  the  foreground  H{0,  31)  represents  collimated 
incident  li^t,  which  is  scattered  into  the  31st  ring  detector,  and 
the  axis  in  the  left  foreground  is  the  sin^e-scattering  signature  for 
a  conventional  difiraction  system  (collimated  input  scattered  over 
31  rings).  The  ridge  at  the  right  of  the  figure  culminating  at 
H(31,  31)  is  caused  by  the  fact  that  the  outer  RSI  ring  detectors  are 
much  larger  than  the  inner  ones.' 
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foreground,  lower  left,  corresponds  to  the  first  col¬ 
umn  of  H,  i.e.,  the  single-scattering  signature  that 
would  be  obtained  with  colUmated  =  0)  incident 
radiation  ^ls  given  by  Eq.  (2).  For  f  =  125  mm  and 
Natl.  Bur.  Stand.  (U.S.)  SRM  1003a,  the  peak  signal 
occurs  in  the  region  of  ring  detectors  10  emd  11.  The 
upper  left  element  H  (0,  0)  corresponds  to  light 
incident  at  @  =  0,  which  is  Uterally  forward  scattered 
so  that  0,tt  =  0.  As  one  moves  through  columns  of  H 
(i.e.,  moving  toward  the  upper  right  background)  the 
matrix  row  at  which  the  scattering  maximum  occurs 
shifts  toward  the  outer  detectors  (larger  i  for  H^). 
Finally  at  H(31,  31)  most  of  the  energy  incident  at  an 
angle  corresponding  to  the  large  31st  ring  of  the  RSI 
detector  is  rescattered  into  the  same  detector.  Plots  of 
H„,  matrices  calculated  made  by  using  Eq.  (7)  for  the 
same  particle-size  distribution  and  optical  system  but 
at  several  optical  depths  are  shown  in  Fig.  4.  The  H 
and  matrices  are  identical  in  the  limit  of  small 
optical  depth  b,  and  the  similarities  are  seen  by 
comparing  Figs.  3  and  4(a).  However,  as  b  increases 
the  first  columns  of  (foreground)  indicate  the 
scattering  signatures  that  would  be  obtained  for 
collimated  incident  radiation  as  given  by  Eq.  (2).  The 
shift  of  the  scattered  energy  toward  larger  angles 
(lower  rows  of  the  matrices)  characteristic  of  multiple 
scattering  is  apparent  in  the  progression  through 
Figs.  4(a)-4(c).  That  this  shift  toward  larger  angles 
resembles  small  particles  to  a  conventional  near- 

(0,31) 


b  =  2.0 


(31.0) 
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(3UI) 


(31.0) 
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forward-scattering  instrument  has  been  well  docu¬ 
mented.’® 

The  inverse  scattering  scheme  proposed  here  re¬ 
quires  measurement  of  the  matrix  H„  and  the  optical 
depth  b  and  the  numerical  determination  of  H  with 
Eq.  (8).  For  this  work  we  first  calculated  H  and  then 
H,„  as  shown  in  Figs.  3  and  4.  A  simulation  of  the 
experimental  measurement  of  H„,  was  then  made  by 
perturbing  the  calculated  elements  of  H,„  with  a 
Gaussian  random  number  generator.  The  random 
noise  added  was  multiphcative  and  was  specified 
whereby  the  standard  deviation  of  the  Gaussian 
random  noise  was  a  fraction  of  the  signal  on  each 
detector. 

Another  experimental  difiiculty  will  be  measure¬ 
ment  of  the  di^onal  elements  of  H„„  which  represent 
light  scattered  into  the  same  detector  aperture  as  the 
incident  radiation.  Since  the  scattering  is  often  an 
order  of  magnitude  or  more  smaller  than  the  unscat¬ 
tered  (transmitted)  energy,  the  elements  will  be 
difficult  to  measure  accurately.  In  the  simulation  of 
the  experiment  the  diagonal  elements  were  deter¬ 
mined  by  averaging  the  ac^acent  off-diagonal  ele¬ 
ments.  Imphcit  in  this  procedure  is  that  optical  and 
electronic  cross  talk  between  detector  apertiu’es  is 
ne^gible. 

After  the  matrix  was  perturbed  to  simulate  an 
experimental  measurement  the  inversion  indicated 
by  Eq.  (8)  was  performed.  The  results  are  presented 


Fig.  4.  Calculated  multiple-scattermg  redistribution  matrix  H„. 
The  redistribution  matrix  H.,  was  calculated  for  the  particle-size 
distribution  of  National  Bureau  of  Standards  (U.S.)  SRM  1003A 
for  a  laser  difiraction  system  that  uses  using  the  RSI  detector 
geometry*’*  with  a  transform  lens  focal  length  of  125  mm  and  the 
successive  order,  discrete  ordinates  multiple-scattering  model.' 
The  plots  correspond  to  three  different  'values  of  the  optical  depth  b 
as  indicated.  The  upper  left  element  of  H„,  i.e.,  H„(0,  0),  is  shown 
on  the  left  and  applies  to  light  incident  at  the  smaUest  angle  6o 
(0,  s:  0  or  collimated  input),  which  is  scattered  back  into  the  same 
an^e  (identically  forward  scattering)  after  passing  throu^  the 
entire  medium.  The  lower  point  in  the  foreground  of  the  figure 
HJ.0,  31)  represents  collimated  incident  light,  which  is  scattered 
into  the  31st  ring  detector,  and  the  axis  in  the  left  foreground 
shows  the  scattering  signature  that  would  be  obtained  by  a 
conventional  diffraction  system  (collimated  input  scattered  over  31 
rings).  The  progression  through  increasing  values  of  b  shows  the 
expected  shift  of  scattered  light  to  larger  angles. 
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Fig.  5.  First  columns  of  the  sin^e-scattering  redistribution  ma¬ 
trix  H  reconstructed  in  numerical  experiments  by  using  an  inver¬ 
sion  of  qmthetic  scattering  data  given  by  Eq.  (8).  The  assumed 
optical  system  had  an  RSI  ring  detector^  with  f  =  125-mm 
transform  lens,  and  the  particle-size  distribution  assumed  was 
National  Bureau  of  Standards  (U.S.)  SRM  1003a.  The  calculated 
elements  of  the  matrix  (shown  in  Fig.  4)  were  perturbed  with 
two  different  levels  of  simulated  Gaussian  (multiplicative)  noise 
with  standard  deviations  of  (a)  1%  and  (b)  5%. 

in  Fig.  5,  which  shows  the  first  column  of  H  as 
obtained  from  Eq.  (8).  To  review,  the  first  column  of 
H  is  an  estimate  of  the  single-scattering  signature, 
which  would  be  obtained  by  using  collimated  on-axis 
incident  radiation.  In  conventional  particle-sizing 


Fig.  6.  RMults  of  the  numerical  inversion  of  the  scattering 
signatures  of  Fig.  5  Reconstructed  particle-size  distributions  are 
shown.  The  size  distributions  were  obtained  by  solving*  Eq.  (9), 
where  the  scattering  signature  8  was  taken  as  the  first  column  of 
the  reconstructed  H  matrices  as  shown  in  Fig.  6. 


methods  using  near-forward  scattering  it  is  the  first 
column  of  H  that  is  hopefully  measured  1^  the 
detectors.  For  that  reason,  then,  inversion  schemes 
used  in  conventional  methods  for  optically  thin  aero¬ 
sols  can  be  used  on  the  first  column  of  the  recon¬ 
structed  H  matrix  to  obtain  the  particle-size  distribu¬ 
tion.  The  various  schemes  for  performing  the  inverse 
single-scattering  problem  are  discussed  elsewhere.®  * 
Tlie  performance  of  the  proposed  method  based  on 
the  numerical  studies  is  demonstrated  in  Fig.  6.  Here 
the  actual  and  reconstructed  size  distributions  for 
SRM  1003a  are  shown  for  several  optical  depths  and 
for  5%  simulated  Gaussian  random  noise.  Reasonable 
reconstructions  of  the  particle-size  distribution  are 
feasible  up  to  6  =  6  for  the  conditions  studied  here. 

Conclusions 

A  scheme  for  solving  the  inverse  scattering  problem 
in  optically  thick  m^a  has  been  presented.  Numeri¬ 
cal  experiments  indicate  that  the  method  is  feasible 
for  reconstructing  particle-size  distributions  in  realis¬ 
tic  measurement  conditions. 
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INTRODUCTION 

The  optical  measurement  technique  commonly  kiwwn  as  the  Fraunhofer  diffraction  particle  sizing  method  involves 
measuring  the  near  forward  scattering  signature  of  an  ensemble  of  large  particles  followed  by  a  mathematica!  inversion 
procedure  to  determine  the  particle  size  distribution  function.  Of  interest  in  this  paper  are  errors  introduced  by  a  lack  of 
knowledge  of  the  particle  relative  refractive  index.  The  near-forward  scattering  signature  of  an  ensemble  of  particles  will 
change  as  the  refractive  index  of  the  particles  approaches  that  of  the  surrounding  medium.  Under  these  conditions  (and 
using  a  geometric  mtics  model)  the  refracted  light  is  increasingly  directed  into  near-forward  angles,  interfering  with  the 
^ffrect^  light,  arxl  thereby  altering  the  near-forward  scattering  pattern.  A  particle  sizing  instrument  based  entirely  on 
Fraunhofer  dif^ction  formulas  will  give  erroneous  results  when  the  refractive  index  of  the  particles  is  similar  to  that  of  the 
surrounding  medium.  It  is  these  changes  in  the  near  forward  scattering  signature  as  the  refractive  indices  change,  and  their 
effects  on  the  size  distribution  which  would  be  measured  by  such  a  particle  sizing  instrument,  which  are  considered  in  this 
paper. 

MATHEMATICAL  FORMULATION 

The  most  general  equation  governing  Fraunhofer  diffraction  particle  sizing  is  a  Fredholm  integral  equation  of  first  order 
and  first  kind  and  can  be  formulated  as: 


1(0)  =  Iinc/k2  /  kb(a.0)  nb(a)  da  (I) 

0 

where:  i(0)  is  the  intensity  (W/sr)  diffracted  at  near-forward  scattering  angles  0;  Ij„(.  is  the  irradiance  (W/m^)  incident  on 
the  particles  (assumed  constant);  k  is  the  wavenumber  2n/7-;  X  is  the  wavelength;  a  is  the  size  parameter  tiD/X  where  D  is 
the  particle  diameter,  n()(a)  s  n(a)afi  where  n(a)da  is  number  of  particles  in  the  laser  beam  with  sizes  between  a  and 
a  +  da  (n  is  an  unnormalized  probability  density);  b  is  a  scaling  parameter  determined  by  the  instrument  designer  which 
represents  the  measure  of  the  quantity  of  particles  used  as  the  solution  basis  (e.g.  for  number,  area,  or  volume  bases  the 
parameter  b  would  equal  0,  2,  or  3  respectively);  and  iq)(a,0)  is  a  general  scattering  function  which  gives  the  scattering 
contribution  of  a  urtit  quantity  of  particles  of  size  a  into  angle  0.  Particle  sizing  using  Eq.  (1)  requires  measurement  of  i(0) 
followed  by  determining  the  nb(a)  which  would  produce  a  calculated  1(0)  signaniie  (using  Eq.  (1))  which  best  fits  in  some 
sense  the  measured  i(0).  Ge^y,  the  function  kt)(a,0)  must  be  known  to  use  Eq.  (1)  in  the  inversion  process,  and  full 
Lorenz-Mie  expressions  or  approximations  such  as  Fraunhofer  diffraction  theory  may  be  used. 

Equation  (1)  represents  an  infinite  set  of  equations  (one  for  each  value  of  0)  in  an  infinite  set  of  unknowns  (nt,(a)  for  all 
possible  a).  In  practice,  however,  it  is  impossible  to  make  an  infinite  number  of  measurements  of  the  scattered  light 
distribution,  and  therefore  a  finite  number  of  measurements  at  a  set  of  discrete  angles  is  usually  made.  By  discretizing  the 
equations  we  obtain  a  system  of  me  equations  in  m^  unknowns  where  me  is  the  number  of  discrete  detectors  and  m^  is  the 
number  of  discrete  size  classes.  We  write  the  linear  system  as: 

I  =  K  •  N  (2) 

In  Eq.  (2)  the  me  elements  of  the  vector  I  represent  the  measured  scattering  signature;  the  m^  elements  of  N  contain  the  bth 
partial  moments  of  the  number  of  particles  in  the  size  class;  and  K  is  the  instrument  matrix  whereby  element  Kjj  represents 
the  diffraction  contribution  of  a  unit  measure  of  particles  in  the  jth  size  class  onto  the  ith  detector.  The  elements  of  K  are: 

Kjj  =  J  kt,(a,0i)  Wo(a,cq)  we(0,0i)  da  d0  (3) 

0 

where  wa(a,aj)  is  a  weighting  function  for  the  jth  size  class  arxl  we(0,0i)  is  a  weighting  function  describing  the  sampling 
process  of  the  deteaor  apertures.  Now  the  solution  or  measured  particle  size  distribution  indicated  by  N  in  Eq.  (2)  can  in 
theory  be  obtained  by  inverting  the  matrix  K. 

When  the  refractive  index  of  the  medium  approaches  that  of  the  particles,  the  assumptions  inherent  in  a  kemal  ba.sed  on 
Fraunhofer  diffraction  are  no  longer  valid.  Calculations  of  near  forward  scattering  signatures  in  this  case  are  based  on 
theory  presented  by  van  de  Hulst  [1957J  where  the  term  anomahus  diffraction  was  coined  to  describe  the  condition  when 
a  is  large  and  the  refractive  itxiex  denoted  m  approaches  that  of  the  surrounding  medium.  A  parameter  p  is  defined  as: 

p  =  2a|m-l|  (4) 

which  physically  represents  the  phase  lag  suffered  ^  a  ray  that  passes  through  a  sphere  along  a  full  diameter  relative  to  a 
ray  whjch  passed  undisturbed  through  the  surrounding  mei^um.  The  presence  of  the  sphere  does  not  alter  the  amplitude  of 
the  field  at  some  position  beyond  the  sphere,  but  rather  affects  only  the  phase.  The  phase  lag  can  be  calculated  at  some 
plane  past  the  sphere  ^  a  direct  application  of  Huygen’s  princqrle  gives  both  the  extinction  and  scattering  diagrams.  The 
intensity  distribution  is  found  by  calculating  the  amplitude  function  A(p.z)  where  z  «  a0.  An  angle  t  is  defined  as  the 
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0 

which  is  derived  by  qiplying  Huygen’s  principle  to  s  plane  just  beyond  the  sphere.  The  following  eiqiansioo  of  the  integal 
anumes  that  m  is  real  and  the  abwipdon  coeffideitt  equals  aro.  The  real  and  imaginary  parts  of  the  integral  A(p^)  are 
evaluated  separately.  If  the  substitution  y  ■  {nfiyx  is  made,  the  imaginary  part  of  the  iiKegtal  is  transfomii^  to  Sonine’s 
second  integ^  and  is  represented  by: 

1mA  ■  (p/y2)  (Ry/2)*^  l3/2(y),  where  y2  «  p2  +  *2  (6) 

The  real  part  is  expreswd  as  aries  erqiansions.  Two  different  expansions  are  necessary  for  small  and  large  values  of  p. 
For  small  p  the  first  aries  expansion  is: 

ReA  -  p2(l/z2)J2(z)  -  (P^/1*3X1/z3)J3(z)  +  (p6/l*3*5Xl/x<)J4(z)  + ...  (7) 

The  expanston  for  large  values  of  p  is  given  by: 

ReA  -  (I/z)J|(z)  +  (p/y2XJty/2)‘/2Y3/2(y)  +  (l/p2)Jo(z)  +  (l*3/p4)zJi(z)  + ...  (8) 

where  Y  is  the  Besal  function  of  the  acond  kind.  The  first  aries  converges  for  any  combination  of  p  and  z  while  the 
acond  converges  only  for  p>z.  Hnally,  the  scattering  function  k|,  in  Eqs.  (1)  aixl  (3)  for  the  arxtmalous  difEraction  caa  is: 

kb(a.e)  -  A2(a0y(ob-2)  (9) 

where  A2  b  (ImAp  +  (ReA)2  as  expresad  by  Eqs.  (6-8).  To  quantify  the  error  of  the  Fraunhofer  diffraction 
rqiproximation  when  anomalous  diffiaction  conditions  are  present  we  tove  generated  near  forward  scattering  sigttatutes 
based  on  anomalous  diffraction  theory  for  MBS  standard  reference  material  (SRM)  l(X)3a. 

RESULTS 

A  plot  of  the  intensity  normalized  over  the  angle  range  0-2.7”  is  shown  in  Fig.  1  for  MBS  l(K>3a,  which  contains  particle 
sias  ranging  from  S-63pm.  As  expected  when  the  value  of  the  refractive  ii^ces  approach  unity,  the  intensity  at  larger 
angles  incream  u  the  refracted  light  interferes  with  the  diffracted  light  in  this  angle  range.  When  the  ratio  of  the 
refraaive  indices  is  greater  than  1.2  the  anomalous  diffraction  results  agree  with  Fraunhofer  diffraction.  It  should  also  be 
noted  that  the  anomalous  diffraction  results  in  this  angle  range  were  verified  by  Lorenz-Mie  theory. 

Rgure  2  is  a  plot  of  the  relative  deteaor  output  vs.  detector  number  for  the  RSI  ring  detector.  Results  are  found  by 
integrating  the  intenrity  distribution  of  Fig.  1  over  each  detector  area.  Refractive  indices  close  to  one  show  more  energy  on 
the  outer  rings.  Again,  recall  that  results  are  for  the  same  particle  size  distribution,  the  only  variable  is  the  refractive  index. 

If  an  instrument  mttering  matrix  K  is  calculated  using  a  Fraunhofer  diffraction  kernel  Iq,  then  sizing  errors  will  result  if 
the  actual  mttering  process  is  not  accurately  modeled  by  diffraction  theory.  The  varying  scattering  signatures  of  Fig.  2 
repre^t  different  1  vectors  in  Eq,  (2),  which  for  the  same  K  must  give  different  measured  size  distributions  repre.sented  by 
N.  Since  small  particles  scatter  light  at  relatively  large  angles,  the  anomalous  diffraction  (refractive  scatter)  contributions 
on  the  outer  detectors  for  small  refractive  indices  in  Fig.  2  appear  to  come  from  small  particles  to  an  instrument  based  on 
I^unbofer  theory.  This  can  be  seen  from  the  reconstruaed  size  distributions  shown  in  Fig.  3  where  a  small  particle  mode 
(in  the  lotrer  sire  classes)  appears  for  the  two  lowest  refractive  indices.  These  resulting  errors  ore  quantified  in  terms  of 
represratative  sire  parameters  as  would  be  measured  an  instrument  based  on  Fraui^fer, diffraction  (Dvx  is  defined 
according  to  ASTM  E799  such  that  the  volume  fraction  of  particles  with  diameters  smaller  than  Dv,i  is  x,  and  avx  » 
n/^Dvx). 

The  data  of  Rg.  4  show  three  regimes,  two  for  the  limiting  cases  of  very  small  and  very  large  relative  refractive  indices  and 
a  third  in  the  intermediate  m  range.  For  very  large  reactive  indices,  the  refracted  light  is  scattered  into  large  mgles 
beyorxi  the  collection  aperture  of  the  ^stem  (see  Rg.  1).  Here  the  Fraunhofer  diffiaction  theory  is  applicable,  and  the 
errors  are  sniall.  For  very  small  refractive  irxlices,  the  scanering  is  very  small  because  the  contrast  between  the  particle  arid 
the  medimn  is  very  weak.  Here,  although  the  refracted  light  does  mix  with  the  near-forward  diffracted  light  the  refracted 
contrilmtirm  is  weiJc  and  again  diffraction  theory  is  reasonable  good  as  demonstrated  by  small  errors.  In  tire  intermediate 
refractive  index  range  (from  about  1.004  to  1.2  in  this  example),  the  refracted  light  contribution  is  significant  arid  falls  imo 
the  collection  aperture.  This  causes  a  dip  in  the  curves  of  Rg.  4,  i.e.  the  refracted  light  entering  Uw  near-forward  angles 
cause  the  insnument  based  on  diffraction  theory  to  undersize  by  as  much  as  70%.  In  summary,  significant  errors  in  particle 
rize  distribution  reconstructions  based  on  Fraunhofer  diffiaction  theory  will  occur  when  the  ratio  of  tlK  refractive  indices  is 
in  an  iiaermediate  range.  Better  approximations  or  the  full  Lorenz-Mie  scattering  theory  should  be  used  in  instruments 
operating  in  these  refractive  index  ranges. 
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